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A  process-based  model  (UFET)  for  deep-submicron  bulk-silicon 
MOSFETs  is  developed  and  verified  with  numerical  device  simulations  and 
measured  data.  The  charge-based  model  is  physical  with  accountings  for  the 
predominant  short-channel  (e.g.,  charge  sharing,  drain-induced  threshold  reduction 
and  velocity  saturation)  and  extremely  scaled-technology  (i.e.,  energy  quantization 
and  polysilicon-gate  depletion)  effects  in  MOSFETs.  The  key  to  UFET  is  the 
characterization  of  the  bias-dependent  two-dimensional  regions  near  the  source/ 
drain  junctions  which  can  extend  over  a  significant  fraction  of  the  metallurgical 
channel  length.  When  these  two-dimensional  regions  near  the  junctions  are  modeled, 
the  physical  charge-sheet  model  can  be  applied  to  the  remaining  "quasi-two- 
dimensional"  channel  length  to  define  the  channel  current  and  terminal  charges, 
without  resorting  to  empiricism  to  account  for  the  short-channel  effects.  Special 


attention  paid  to  continuity  in  the  derivation  of  the  model  formalism  yields  a  physical 
C^  model  applicable  to  analog  and  digital  CMOS  circuit  design.  The  small  number 
of  physical,  process-based  parameters  simplifies  the  mode!  calibration,  and  renders 
the  model  suitable  for  predictive  device/circuit  simulation,  statistical  simulations 
and  circuit  sensitivity  analyses  based  on  known  or  presumed  process  variations. 
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CHAPTER  1 
INTRODUCTION 


Bulk-silicon  CMOS  technology  continues  to  be  predominant  in  VLSI 
applications  due  largely  to  the  prohibitive  capital  equipment  costs  required  for  any 
large  shift  in  processing  materials  or  procedure.  Although  it  has  been  challenged  by 
GaAs  with  its  higher  carrier  mobility  and  silicon-on-insulator  (SOI)  with  its  reduced 
junction  capacitances  for  more  than  a  decade,  it  has  defied  predictions  that  these  new 
technologies  would  soon  dominate;  the  reality  is  that  these  two  new  technologies 
remain  confined  to  niche  markets  at  this  time. 

With  this  in  mind  it  is  important  to  develop  accurate  compact  models  for 
extremely  scaled  bulk  CMOS  devices  for  use  in  process  development/integration  and 
statistical  modeling  as  well  as  circuit  design.  Technology  CAD  (TOAD)  requires 
integrated,  physics-based,  CPU-efficient  device  models  for  predictive  circuit 
simulation.  Although  many  useful  MOSFET  models  are  available,  they  generally 
lack  an  adequate  physical  basis  and/or  acceptable  computational  efficiency,  in 
addition  to  containing  discontinuities  in  their  formulations.  With  a  conventional 
MOSFET  model  such  as  BSIM  [Che96],  a  large  number  of  empirical  device  model 
parameters  are  optimized  based  on  measured  data  (e.g.,  via  least-squares  fitting)  for 
circuit  simulation.  The  needed  optimization  misses  correlations  between  device 
model  and  technology  parameters  which  may  yield  non-unique  (erroneous)  circuit 
simulations. 

1 


However,  the  use  of  application-specific  analytic  TCAD  tools  has  been 
suggested  as  a  viable  way  of  attaining  the  desired  computational  efficiency  [Mar88] 
while  maintaining  a  strong  physical  basis.  This  approach  is  exploited  through  the 
incorporation  of  application-specified  (i.e.,  a  predetermined  device  structure) 
analytic,  but  semi-numerical,  device  models  in  the  circuit  simulator.  Unlike 
conventional  optimization-based  parameter  extraction,  the  evaluation  of  the  model 
parameters  is  based  on  an  analysis  of  the  device  structure  using  criteria  consistent 
with  the  device  physics.  Therefore,  the  parametric  correlations  are  inherently 
accounted  for,  and  the  sensitivities  of  device  and  circuit  performance  to  variations  in 
the  fabrication  process  can  be  reliably  and  efficiently  predicted. 

A  bipolar  junction  transistor  (BJT)  model  that  meets  the  criteria  required 
for  pragmatic  semi-numerical  device/circuit  simulation  was  developed  previously 
[Jeo89,  Hon94  and  Cha97].  The  model  is  based  on  linked  regional  analyses  of  the 
device  structure,  and  provides  a  physical  alternative  to  the  widely  used  empirical 
Gummel-Poon  model  [Get78].  With  its  implementation  in  SPICE2  [Nag75]  and 
incorporation  of  the  parameter  evaluator  SUMM  [Gre90],  a  viable  mixed-mode  BJT 
device/circuit  simulator  was  created. 

The  deficiencies  of  existing  compact  short-channel  MOSFET  models 
preclude  predictive  simulations  and  use  of  the  models  in  this  pragmatic  approach.  As 
MOSFET  channel  lengths  have  been  scaled  into  the  deep-submicron  region, 
conventional  short-channel  circuit  models  have  generally  been  developed  via 
empirical  extensions  of  physical  long-channel  models.  These  extensions  were  based 
on  empirical  modifications  to  simple  one-dimensional  characterizations  of  the 


electric  field  to  account  for  the  short-channel  effects,  which  are  controlled  by 
complex  two-dimensional  field  distributions.  Most  MOSFET  models  do  not 
explicitly  account  for  the  various  short-channel  effects;  for  example,  BSIM  relies 
heavily  on  empirical  "fitting"  expressions  that  are  intended  to  account  for  these 
effects  without  distinguishing  them.  Such  approaches  were  taken  to  preserve  the 
accuracy  and  efficiency  required  for  simulation  of  large  circuits,  but  with  the  loss  of 
predictive  capability.  In  order  to  achieve  predictive  simulation,  rigorous  two- 
dimensional  analyses  of  the  electric  field  are  necessary,  but  the  required  computation 
expense  is  undesirable  for  circuit  simulation.  This  dissertation  presents  a  process- 
based  model  for  deep-submicron  MOSFETs  (UFET)  that  relies  on  approximate  two- 
dimensional  analyses,  which  are  computationally  efficient,  to  account  for  short- 
channel  effects,  viz.,  threshold  voltage  reduction  due  to  charge  sharing  [SzeSl], 
drain-induced  charge  enhancement  [Vee88],  drain-induced  barrier  lowering  [Tro79] 
and  channel-length  modulation  [SzeSl].  These  solutions  of  two-dimensional 
analyses  are  physically  linked  and  are  expressed  in  simple  terms  that  yield  physical 
insight  into  device  behavior  as  well  as  a  predictive  circuit  model  with  utility  for 
TCAD,  e.g.,  for  statistical  simulations  and  circuit  sensitivity  analyses,  as  well  as  for 
design  in  extremely  scaled  technologies. 

Chapter  2  addresses  extremely  scaled  MOSFET  modeling  issues,  and 
presents  an  overview  of  the  previous  generation  of  bulk  MOSFET  models 
highlighting  the  problems  associated  with  regional  analyses.  A  refined  iteration 
scheme  for  properly  modeling  LDD  devices  is  also  documented. 


Chapter  3  presents  the  new  process-based  UFET  circuit  model  which 
includes  quasi-two-dimensional  analyses  that  do  not  have  discontinuous  regional 
boundaries.  The  model  applies  directly  to  extremely  scaled  MOSFETs  with 
retrograde  channel  and  extended  source  and  drain  regions,  and  comprehensively 
accounts  for  short-channel  effects,  field-dependent  mobility,  impact  ionization  and 
device  parasitics.  Iterative  solution  techniques  are  required  in  the  model  formalism 
to  predict  channel  current,  impact-ionization  current  and  quasi-static  regional 
charges  (to  account  for  transient  charging  currents)  from  technology  parameters  with 
minimal  reliance  on  empirical  modeling. 

Chapter  4  overviews  the  intrinsic  model  algorithm  and  the  parasitic/ 
extrinsic  portion  of  UFET.  It  also  details  the  implementation  of  the  new  model  (in 
both  MMSPICE  and  a  platform-independent  Application  Program  Interface)  and  a 
unique  procedure  for  process-based  parameter  evaluation  which  focuses  on  a  limited 
number  of  parameters  in  each  key  region  of  operation.  This  chapter  presents 
verification  of  the  model  based  on  device  simulations  and  measured  data  from 
0.15|J,m  and  0.25|Xm  CMOS  technologies. 

Chapter  5  presents  a  simulation  exercise  in  a  0.1  |im  technology  with 
MEDICI  which  highlights  the  ease  of  parameter  evaluation  for  this  model  and  the 
power  of  a  process-based  model  with  strong  correlation  to  process  variations.  Also 
presented  are  the  results  of  a  sensitivity  study  using  the  mentioned  0.1 5nm 
technology  for  the  fabrication  of  a  two-stage  folded  cascode  operational  amplifier. 

Chapter  6  summarizes  this  work  and  presents  suggestions  for  model 
extensions  and  future  research  in  this  area. 


CHAPTER  2 
SCALED  MOSFET  MODELING  ISSUES 


2.1   Introduction 

As  CMOS  technologies  are  being  scaled  to  sub-O.Ifi,m  gate  lengths,  it  is 
becoming  apparent  that  conventional  methods  of  MOSFET  modeling  for  circuit  simulation 
are  inadequate.  For  example,  the  need  for  model  continuity  [Haj79]  in  all  regions  of 
operation  is  inferred  when  designing  circuits  for  low-voltage  digital  and  analog  operation, 
which  accompanies  extreme  scaling.  Further,  lateral  as  well  as  vertical  variations  in 
channel  doping  for  controlling  (two-  and  three-dimensional)  short-channel  effects  [  Yan9 1  ] 
imply  the  need  for  truly  physical  modeling,  which  should  not  rely  on  global  optimization 
of  model  parameters.  BSIM3  [Che96]  is  currently  the  most  widely  used  model  for  circuit 
simulation,  but  its  empirical  nature,  as  reflected  by  the  large  number  of  parameters  that 
require  copious  amounts  of  data  for  parameter  extraction,  portends  limits  to  its  utility  for 
extremely  scaled  technologies.  Although  an  empirical  MOSFET  model  can  be  used  for 
circuit  design  in  a  mature  CMOS  technology,  it  is  necessary  to  have  a  physical  model 
to  facilitate  the  design  of  new  semiconductor  processes  and  the  modeling  of  how  the 
devices  behave  under  normal  process  variations,  as  well  as  to  ease,  or  even  enable 
the  parameter-extraction  process.  Further,  such  a  physical  model  may  be  essential 
even  for  circuit  design  in  mature  extremely  scaled  technologies. 

In  this  chapter,  a  physical  model  for  the  bulk  MOSFET  (QBMOS) 
implemented  in  MMSPICE-3.0  [Fos94]  will  first  be  reviewed.  Then  results  from  a 


technology  transfer  with  AT&T  will  be  discussed  in  terms  of  general  implications  for 
all  MOSFET  modeling,  including  in  SOI.  A  discussion  of  model  refinements  leading 
towards  applications  in  extremely  scaled  technologies  is  also  presented  as  a  lead-in 
to  Chapter  3. 

2.2  QBMOS:  A  Physical  Bulk  MOSFET  Model 
QBMOS  [Gre93a]  is  a  significant  extension  of  the  previously  developed 
short-channel,  charge-based  model  for  thin-film  SOI  MOSFETs  [Vee88].  The 
noteworthy  extensions  are  accountings  for  lightly  doped  drain/source  (LDD/LDS) 
and  impact  ionization  [Cho91]  as  well  as  subthreshold  current  [Gre93b]  in  the  device 
structure  shown  in  Figure  2.1.  The  new  four-terminal  (source,  drain,  gate  and 
substrate)  model  is  based  on  the  assumptions  of  uniform  channel  doping  with  a 
constant  depletion  depth  along  the  channel.  The  analysis  of  the  carrier-transport 
mechanisms  in  the  device  requires  an  iterative  numerical  solution  mainly  due  to  the 
LDD  analysis,  which  must  be  physically  linked  to  the  channel  analysis.  The  resulting 
semi-numerical  model  algorithm  is  implemented  in  SPICE2  source  code  creating 
MMSPICE-3.0,  enabling  efficient  mixed-mode  device/circuit  simulation  since  the 
key  model  parameters  are  either  structural  or  physics-based. 

QBMOS  is  based  on  a  regional  approach  with  physical  modeling  of  the 
channel  current  in  weak  (dominant  diffusion  current)  and  strong  (dominant  drift 
current)  inversion.  This  regional  approach  simplifies  the  channel-current  modeling 
in  the  two  regions,  but  requires  a  cubic  spline  to  join  these  analyses  across  the 
moderate-inversion  region,  which  has  physically  boundaries  (Vjw  and  V-j-j)  [Tsi82]. 
In  devices  where  short-channel  effects  significantly  affect  the  output  characteristics 
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Figure  2. 1   Cross-sectional  view  of  a  modern  bulk  MOSFET  with  LDD  and  LDS. 


of  the  device,  it  is  necessary  to  account  for  the  two-dimensional  effects  on  surface 
potential  and  inversion  charge.  The  accountings  for  charge-sharing  and  drain- 
induced  effects,  with  the  LDD  option,  are  now  developed  to  complete  the 
presentation  of  the  model. 

2.2.1   Charge-Sharing 

Short-channel  devices  exhibit  two-dimensional  effects  for  Vpj  =  0  due 

to  the  lateral  electric  fields  associated  with  the  source  and  drain  junctions.  This  effect 
is  handled  differently  in  strong  and  weak  inversion. 

In  strong  inversion,  a  geometrical  charge  sharing  model  [Yau74]  is  used 
to  define  an  effective  gate-controlled  body  charge  (areal)  density 


H  /,..'^_   ' 


'-rlJ'^^1?-',, 


(2.1) 


Vb(eff)   - 

■■?',-.■  ' 

where  Xj  is  the  source/drain  junction  depth  and  x^^  is  the  one-dimensional  depletion 
depth.  Reducing  the  effective  body  charge  causes  the  gate-controlled  inversion 
charge  density  to  increase  through  Gauss's  law  (viz.,  the  threshold  voltage  is 
reduced)  with  the  accompanying  increase  in  channel  current. 

The  weak-inversion  effect  of  charge  sharing  is  modeled  more  physically 
by  solving  a  differential  equation  resulting  from  Gauss's  law  in  the  regions  near  the 
source  and  drain  [Gre93b]: 


c   +i^ 

dy  EsXj„  2 


By  modeling  the  length  of  these  regions,  a  gradual  channel  can  be  defined  where  the 
dominant  diffusion  current  can  be  modeled  using  one-dimensional  theory.  A  self- 
consistent  assumption  for  the  incremental  field  at  the  boundaries  of  these  regions 

would  be   — —  »  AE     where   0p„   is  the  electron  quasi-Fermi  potential,  but  for 

simplicity  this  boundary  condition  is  approximated  as  AE  (L  -  L^)  h  0 .  The  resulting 
analysis  gives  the  length  of  the  gradual  channel 


(2.3) 


to  be  used  to  evaluate  the  subthreshold  channel  current  with  an  accounting  for  charge 
sharing. 

2.2.2  Drain-Induced  Effects 

In  short-channel  devices,  applied  drain  biases  enhance  the  channel  current 
through  Drain-Induced  Barrier  Lowering  (DIBL)  [Tro78]  in  weak  inversion  and 
Drain-Induced  Charge  Enhancement  (DICE)  [Vee88]  in  strong  inversion.  DIBL 
results  from  the  direct  modulation  of  the  surface  potential  (\)/j,)  by  Vq^,  and  DICE  is 
the  indirect  modulation  of  the  inversion  charge  implied  by  the  two-dimensional 
Poisson's  equation  under  the  gate. 
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Application    of   Poisson's   equation    in    the   depletion    region    can    be 
simplified  to  Laplace's  equation  for  the  incremental  change  in  potential  due  to  Vqs: 


V  d^^d^^o_  (2.4) 

dx^        dx^ 


which  can  be  solved  by  a  quasi-separation  of  variables: 


d!M,_d!^s_„,  (2.5) 

dx-'  dx-^ 


2Vns 
where  r\  =  — f^  [Vee88].  This  first-order  extrapolation  reduces  to  the  gradual- 
channel  approximation  [SzeSl]  for  long-channel  devices. 
Integrating  (2.5)  over  x  yields 


.„  ,  ,       A\|/,(y)-Av(Xj„,y)     r\x^„ 

AE,(y)  = —.  (2.6) 


This  expresses  the  change  in  the  surface  transverse  electric  field  as  a  function  of  the 
change  in  the  potentials  and  the  drain  bias.  Assuming  vertical  fields  in  the  oxide, 
Gauss's  law  is  applied  at  the  Si  -  SiOj  interface  with  (2.6)  to  get 

AQc(y)  =  (c„,  +  ^y^Vsiy)  -  ^-^  ■  (2.7) 
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In  weak  inversion,  it  can  be  shown  that  the  inversion  charge  term  is  much 
smaller  than  the  terms  on  the  right-hand  side  of  (2.7),  so  the  effect  the  drain  bias  is 
seen  as  a  y-independent  increase  in  the  surface  potential  (DIBL): 


AypiBL  =  _!2!l^^do_  ^2.8) 

2{-  '-^ 


<^-3 


The  corresponding  effect  in  strong  inversion  is  seen  as  a  y-independent  increase  in 
the  inversion  charge  along  the  entire  channel  (DICE): 


AQDICE  =  _!i!]^.  (2.9) 


2.2.3  LDD  Analysis 

The  voltage  drop  across  the  LDD  (Vloq)  is  computed  from  the  field 
distribution  in  the  region  [Cho91].  Depending  on  the  maximum  electric  field,  there 
are  three  possible  conditions  in  the  LDD  as  shown  in  Figure  2.2.  Neutral  regions  are 
characterized  by  constant  fields  and  represent  simple  ohmic  drops.  Depletion  regions 
induced  in  the  LDD  by  the  high  channel  field  are  characterized  by  a  constant  slope 
assuming  the  electron  charge  is  negligible, 

-L^  =  -      ^""  =  -S.  (2.10) 

dy  e, 

where  Nldd  '^  ^^^  doping  density  in  the  LDD. 
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Gate 


Drain 


Figure  2.2  Possible  lateral  electric  field  distribution  in  an  LDD  MOSFET  as  a 
function  of  Vps  [Cho91]. 
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The  lateral  electric  field  at  the  body-LDD  junction  is  calculated  using  the 
channel-length  modulation  analyses,  but  it  can  easily  be  shown  [Vee88]  that  it  can 
be  accurately  approximated  as 


o          ^DS  ~  ^LDD  ""  ^DS(eff)  /tun 

E„S j '—^  (2.11) 


where  Vps^^ffj  is  a  result  of  the  channel  analysis  and  I;,  is  the  characteristic  length 
resulting  from  the  analysis  of  the  two-dimensional  velocity  saturated  region. 

The  calculation  of  Vldd  '^  done  as  follows.  For  ( Vp^  -  V^^pp)  <  Vp^.^ff- , 
the  maximum  electric  field  is  assumed  to  be  small  so  that  only  the  ohmic  drop  in  the 
LDD  needs  to  be  considered,  which  corresponds  to  field  distribution  '  I '  in 
Figure  2.2.  Therefore, 

^LDD  =  'ch'^ldd-  (2.12) 

^^  (^ds~^ldd)>  ^DS(eff)'  ^"  """^^  f'^l'^  distributions  '1',  '2'  and  '3'  are  possible 
depending  on  the  value  of  the  maximum  field.  If  E^  <  — ,  then  case  '1' 


represents    the    field    distribution    and    (2.12)    is    used    to    compute    Vlod-    If 

CH     LDD 

0<E|^ <SLlpp,  condition  '2'  exists,  and  the  voltage  drop  across  the 


LDD  is  calculated  using 
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2S 


Pm 1 I 

'""      •  +  IchRldd-  (2.13) 


For  conditions  when  E^ >SLLpp,  the  entire  LDD  is  depleted  and  case 

^LDD 

'3'  is  applicable  with  Vldd  given  as 

VLDD   =   LLDD(Em-^^)-  (2.14) 

An  iterative  approach  is  necessitated  by  the  coupling  between  the  channel 
analysis  and  the  voltage  drop  across  the  LDD  through  the  maximum  lateral  electric 
field  at  the  body-LDD  junction  (E„).  In  QBMOS,  the  LDD  drop  is  solved  using  a 

secant  method  with  initial  iterates  of  V^pp  =  0  and  Vj^p  =  0.1  Vp^. 

2.2.4  Impact  Ionization 

QBMOS  accounts  for  generation  current  due  to  (weak)  impact  ionization 
involving  electrons  traversing  the  high-field  region  near  the  drain: 

lGi  =  (M-l)IcH  (2.15) 

where  a  local-field  model  gives  the  ionization  rate  as 

J' 
(M-1)  =  Ja„e  ®dy  (2.16) 

with  a„  and  p^  being  constants. 
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In  the  channel,  the  multiplication  factor  in  (2.16)  can  be  evaluated  as 


[Vee88] 


(M-l)cH  =  p^[VDS-VDS(eff)]e"--^-'"".  (2.17) 

The  generation  current  in  the  LDD  is  evaluated  by  numerical  integration  of  the  linear 
field  distribution  in  the  depleted  LDD  region  [Cho91]: 

a        ^ 
(M-l)LDDSyjel'=ld|E|.  (2.18) 

The  sum  of  (2.17)  and  (2.18)  is  the  total  multiplication  factor  used  to  compute  Iqj  in 
(2.15). 

2.3  Issues  Regarding  Model  Utility 

In  November  1993,  a  two-week  technology  transfer  took  place  between 
the  University  of  Florida  and  AT&T  (Allentown,  Pennsylvania).  This  consisted  of 
installing  MMSPICE-3.0  and  applying  QBMOS  to  a  0.9|im  CMOS  process  as  well 
as  running  a  series  of  thirteen  general  model  tests  similar  to  those  presented  by 
Tsividis  and  Suyama  [Tsi94].  These  tests  are  designed  to  highlight  flaws  in  model 
physics,  modeling  approaches  and  model  implementations  by  qualitatively 
examining  the  output  characteristics  under  different  DC  and  AC  bias  conditions. 

The  necessity  of  model  continuity  generally  precludes  using  regional 
analyses  due  to  problems  introduced  by  the  use  of  spline  functions.  This  problem 
with  the  regional  approach  is  apparent  in  the  QBMOS  transconductance  plot  of 
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Figure  2.3  Example  of  cubic  spline  problem. 
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Figure  2.3.  Since  the  cubic  spline  is  based  on  the  value  of  the  channel  current  and  the 
transconductance  at  the  boundaries,  the  value  of  the  channel  current 
(transconductance)  across  moderate  inversion  does  not  necessarily  reflect  the 
physical  current  (transconductance)  but  only  reflects  the  function  which  is  necessary 
to  satisfy  the  boundary  conditions.  The  cubic  spline  also  introduces  discontinuities 
in  the  higher-order  derivatives  of  the  channel  current  at  the  boundaries  since  the 
cubic  spline  only  guarantees  the  continuity  of  the  function  and  its  first-derivative. 
This  is  apparent  in  the  transconductance  plot  of  Figure  2.4  where  the 
transconductance  is  continuous  but  the  higher-order  derivatives  of  the  channel 
current  are  obviously  not.  A  similar  discontinuity  is  also  common  at  the  onset  of 
velocity  saturation  due  to  the  difference  in  the  channel  current  models  above  and 
below  saturation.  This  shows  up  as  a  discontinuity  in  the  QBMOS  output 
conductance  as  shown  in  Figure  2.5.  The  use  of  a  regional  approach  is  even  further 
strained  when  effects  (e.g.,  charge  sharing)  are  not  modeled  with  the  same  approach 
in  weak  and  strong  inversion.  In  QBMOS,  this  introduces  a  problem  for  the  moderate 
inversion  case  since  the  weak-moderate  inversion  boundary  is  not  a  function  of  the 
charge-sharing  in  the  device  but  the  moderate-strong  inversion  boundary  is.  In  the 
extreme  case,  the  moderate-inversion  region  collapses  to  zero  width  and  the  channel 
current  is  then  discontinuous  at  this  transition  as  shown  in  Figure  2.6. 

Unlike  the  moderate-inversion  channel  current,  the  moderate-inversion 
terminal  charge  modeling  in  QBMOS  is  done  with  a  linear  interpolation  of  the 
terminal  charges  at  the  weak-  and  strong-inversion  boundaries.  While  the  integrated 
charge  across  moderate  inversion  is  correct  with  this  approach  and  adequate  for 
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Figure  2.4  Example  of  transconductance-discontinuity  problem. 
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Figure  2.5  Example  of  output  conductance-discontinuity  problem. 
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Figure  2.6  Example  of  charge-sharing  problem. 
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digital  simulations,  the  moderate-inversion  capacitance  is  a  constant  due  to  the  use 
of  linear  interpolation  which  is  problematic  for  analog  applications.  The  constant 
capacitance  problem  is  exemplified  in  Figure  2.7  which  shows  Cqq.  Also,  note  the 
problem  of  a  constant  capacitance  in  strong  inversion  for  V^s  =  OV  due  to  the 
assumption  of  a  pinned  surface  potential  in  strong  inversion. 

The  regional  approach  has  several  deficiencies  with  regard  to  the 
continuity  of  the  output  characteristics  with  respect  to  the  terminal  voltages. 
Although  discontinuities  such  as  those  found  in  QBMOS  are  not  necessarily  a 
problem  for  digital  simulations,  the  necessity  of  continuous  derivatives  for  analog 
applications  does  not  facilitate  the  regional  approach  for  use  in  extremely  scaled 
technologies.  The  discontinuities  also  have  an  adverse  effect  on  the  overall 
convergence  and  convergence  speed  of  the  model.  The  regional  approach  can  also 
lead  to  additional  problems  when  accountings  for  the  same  effect  differ  significantly 
in  strong  and  weak  inversion. 

Although  the  present  version  of  QBMOS  is  adequate  for  digital 
simulation,  it  is  desirable  to  have  a  model  which  is  applicable  for  both  digital  and 
analog  applications.  To  accomplish  this,  QBMOS  needs  to  be  improved  by  removing 
all  discontinuities  in  the  current  and  terminal  charge  models  to  result  in  a  C^ 

continuity  model  (i.e.,  the  currents,  charges  and  all  derivatives  are  continuous). 
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Figure  2.7  Moderate  inversion  capacitance  problem. 
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2.4  Algorithm  Refinements 

2.4.1  Improved  LDD  Loop 

A  significant  time  penalty  for  the  intrinsic  model  arises  from  the  iterative 
numerical  solution  required  by  the  presence  of  a  non-ohmic  voltage  drop  across  the 
LDD  region.  In  QBMOS,  the  LDD  drop  is  iteratively  solved  using  a  secant  method 
with  initial  iterates  of  V/^^q  =  0  and  V^pp  =  0.1  Vpj  for  the  bias  point  solution. 
Although  this  method  converges  to  the  correct  solution,  the  lack  of  more  precise 
initial  iterates  increases  the  CPU  time  need  to  find  the  solution. 

However  with  (2.11),  the  analysis  can  be  simplified  by  substituting  this 
approximation  for  the  maximum  electric  field  into  (2.13),  yielding 


r^DS  ~  ^LDD  ~  ^DS(eff)  _  [cH^LDD~|^ 

+  IchRldd-  (2.19) 


^LDD     -I 


LUU  2S 

Now,  (2.19)  has  V^dd  °"  ''°'''  sides,  and  through  algebraic  manipulation  it  can  be 
rearranged  to  ^       -  \ 


VgpD  -  4^05  -  VDS(eff)  -  '^"    "-"^  '  -  S1?1Vldd  •  (2-20) 

'  +  2SIchRldd1"?  =  0 


'LDD 

"^DS~^DS(eff)      'chRlDdI^ 
'c  '.'LDD     J 


which  can  be  easily  solved  to  find  the  next  iterate.  This  solution  method  exploits  the 
fact  that  no  quantity  in  (2.19)  is  a  strong  function  of  Vldd,  which  is  a  good 
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approximation  across  variations  in  channel  length  and  bias  conditions.  It  is  also 
possible  to  apply  the  same  approach  to  (2.14),  which  results  in 


rVn<;  -  V,  Dp- VDS(eff)  _  ^^L 

I  2 


V            _i           I   'DS~  "LDP"  ''DS(eff)      ■^^LDD  /tt,n 

'LDD   -   1-LDD   i ;; I-  (^-^i) 


Since  this  equation  is  linear,  it  is  easy  to  solve  for  the  next  iterate,  which  is 


L, 


rVps  ^DS(eff)    ^^lddI 


•c  2       J 

Vldd  =  1 •  (2.22) 


is  done,  resulting  in  values  for  Iqjj  and  Vos^^ffj  which  lead  to  a  value  for  Vldd 
through  either  (2.12),  (2.13)  or  (2.14)  as  determined  by  the  depletion  condition  in  the 
LDD  region.  If  the  LDD  drop  consists  of  a  non-ohmic  component,  the  next  iterate  for 
^LDD  is  then  determined  using  either  (2.20)  or  (2.22),  again  based  on  the  depletion 
condition  in  the  LDD.  This  fixed-point  iteration  scheme  generally  converges  in  two 
or  three  iterations,  offering  a  significant  improvement  over  the  previously  used 
secant  method. 

2.4.2  Storage  of  Initial  Iterates 

Due  to  the  implicit  solutions  comprising  QBMOS,  the  terminal 
characteristics  do  not  have  explicit  partial  derivatives  (which  are  needed  to  form  the 
Jacobian  matrix),  and  this  requires  the  use  of  finite  difference  approximations  to 
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define  the  partial  derivatives.  This  leads  to  the  intrinsic  model  routine  being  called 
four  times,  once  for  the  bias  point  and  three  times  for  the  derivatives  with  respect  to 
gate,  substrate  and  drain  biases.  In  QBMOS,  the  intrinsic  model  calculates  the 
solution  of  the  iterative  Vldo  analysis  using  the  same  initial  iterate  for  each  call. 
However,  the  bias  point  solution  for  the  voltage  across  the  LDD  can  be  stored  and 
used  as  the  initial  iterate  for  the  subsequent  finite  difference  calls  to  the  intrinsic 
model  routine.  The  use  of  the  better  initial  iterate  allows  the  solutions  for  the 
perturbed  biases  to  be  found  in  fewer  iterations  decreasing  the  cpu  time. 

The  use  of  this  initial  iterate,  coupled  with  the  implementation  of  the 
refined  LDD  loop,  yields  a  large  speed  improvement  for  simulations  of  devices  with 
significant  LDD  drops.  As  shown  in  Figure  2.8,  the  original  version  of  QBMOS  takes 
2.3  times  as  long  to  simulate  as  compared  to  a  version  of  model  that  incorporates  the 
previously  discussed  enhancements.  The  exemplary  simulation  was  a  DC  drain 
characteristic  over  6  Vq^'s  (0.5,  I,  2,  3,  4  and  5V)  using  6  different  channel  lengths 
(0.5,  0.7,  0.9,  2,  5  and  10  |im).  The  simulation  included  a  total  of  50  V^s  bias  points 
for  each  Vqs  and  channel  length  for  a  total  of  1800  bias  points. 

2.5  Summary 
MOSFET  models  based  on  regional  analyses,  which  rely  on  a  spline 
functions  to  join  these  regions,  tend  to  have  problems  with  function  derivatives  in  the 
transition  regions  and  also  discontinuities  of  the  higher-order  derivatives  at  the 
boundaries  of  these  regions.  These  problems  imply  the  need  for  MOSFET  models 
which  are  not  based  on  regional  analyses  and  have  C^  continuity.  The  continuity 
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issues  should  be  addressed  in  the  derivation  of  the  model  formalism  to  guarantee  that 
empirical  smoothing  functions  are  not  used  to  produce  the  desired  terminal 
characteristics  as  done  in  many  MOSFET  models.  Then  the  physical  integrity  of  the 
model,  with  C^  continuity,  can  be  preserved. 

As  CMOS  technologies  are  being  scaled  to  sub-0.1|im  gate  lengths,  bulk- 
silicon  devices  are  moving  away  from  uniformly  doped  substrates  to  nonuniform 
vertical  and  lateral  doping  profiles,  but  many  MOSFET  models  are  based  on  a 
uniformly  doped  substrate.  MOSFET  models  for  extremely  scaled  technologies  must 
account  for  the  variations  in  both  the  vertical  and  lateral  profiles,  and  the  effects 
associated  with  the  large  vertical  fields  (i.e.,  polysilicon-gate  depletion  and  energy 
quantization)  found  in  these  devices  while  maintaining  C^  continuity. 

A  new  MOSFET  model  is  needed  which  accounts  for  the  physical  effects 
and  nonuniform  doping  profiles  found  in  extremely  scaled  technologies.  The  model 
should  be  developed  avoiding  regional  analyses  due  to  the  inherent  continuity 
problems  associated  with  this  approach  as  shown  in  this  chapter.  It  should  also 
incorporate  the  improvements  made  to  QBMOS,  i.e.,  the  improved  LDD  iteration. 
The  resulting  model  must  reliably  simulate  deep-submicron  devices,  and  enable 
predictive  circuit  simulation  as  noted  in  Chapter  1. 
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CHAPTER  3 

UFET:  A  PROCESS-BASED  MOSFET  MODEL  WITH  PHYSICAL 

ACCOUNTINGS  FOR  SHORT-CHANNEL  EFFECTS 


3.1  Introduction 

Almost  twenty  years  ago,  the  charge-sheet  model  was  developed  [Bre78] 
for  (long-channel)  MOSFETs  as  a  physical  but  more  compact  alternative  to  the 
model  of  Pao  and  Sah  [Pao66],  which  requires  numerical  solutions  of  integral 
equations.  The  charge-sheet  model,  applicable  in  all  regions  of  operation,  is  based  on 
the  gradual-channel  approximation  [SzeSl]  and  is  accurate  for  long-channel  devices 
where  short-channel  effects  are  negligible  and  the  saturation  of  the  drain  current  with 
increasing  drain  voltage  can  be  modeled  with  the  pinch-off  approximation.  However 
for  advanced  scaled  devices,  short-channel  effects  cannot  be  ignored,  velocity 
saturation  is  important  and  the  lengths  of  the  two-dimensional  regions  near  the 
source/drain  junctions  can  be  a  significant  fraction  of  the  metallurgical  channel 
length.  These  effects  cannot  be  unequivocally  accounted  for  in  the  classical  charge- 
sheet  formalism,  and  hence  it  is  not  used  much  today.  However,  if  the  two- 
dimensional  regions  near  the  junctions  are  modeled  well,  the  physical  charge-sheet 
model  can  be  applied  to  the  remaining  channel  length  to  define  the  channel  current 
without  having  to  resort  to  empiricism  to  account  for  the  short-channel  effects.  The 
development  of  a  new  scaled  MOSFET  model  (UFET)  based  on  this  concept  with 
accountings  for  the  effects  of  threshold  voltage  reduction  due  to  charge  sharing 


28 


29 

[SzeSl],  drain-induced  conductivity  enhancement  [Vee88],  drain-induced  barrier 
lowering  [Tro79],  carrier-velocity  saturation  and  channel-length  modulation  [SzeSl] 
is  described  in  this  chapter. 

Figure  3.1  shows  the  typical  structure  of  the  deep-submicron  enhancement 
mode  (n-channel)  bulk-silicon  MOSFET,  with  lightly  doped  drain  (LDD)  and  source 
(LDS)  regions  and  retrograded  channel  doping.  Such  a  regional  structure  is  designed 
to  control  most  short-channel  effects.  The  original  purpose  of  the  LDD  was  to 
improve  device  reliability  by  reducing  peak  electric  fields  in  technologies  designed 
to  work  with  large  supply  voltages.  However  in  advanced  devices  designed  for  use 
with  reduced  power  supplies,  these  regions  are  more  appropriately  referred  to  as 
shallow  source/drain  extensions,  having  doping  concentrations  significantly  higher 
than  those  that  yield  significant  reductions  in  peak  fields,  and  are  used  mainly  for 
added  control  of  short-channel  effects.  The  presence  of  the  p"*^  region  in  the 
retrograded  structure  reduces  the  short-channel  effects  [Yan9l]  as  well  by  limiting 
the  depletion  depth.  The  new  model  for  this  device  is  simplified  by  representing  the 
retrograde  by  a  low-high  junction,  and  then  characterizing  the  intrinsic  and  extrinsic 
regions  separately,  incorporating  the  LDD  analysis  by  a  fixed-point  iterative  solution 
technique  as  discussed  in  Chapter  2. 

In  the  new  UFET  model,  the  channel  current  in  the  intrinsic  device  is 
characterized  based  on  modeling  both  drift  and  diffusion  for  all  bias  conditions, 
which  eliminates  the  need  for  the  moderate-inversion  spline  function  and  enhances 
the  continuity  of  the  model.  The  channel  current  is  characterized  by  integrating  the 
drift  and  diffusion  components  over  a  region   in  the  channel   where  the  two- 
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Figure  3. 1  The  retrograde  n-channel  MOSFET  structure  with  LDD/LDS  regions  and 
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dimensional  effects  are  not  large,  but  are  merely  "quasi-two-dimensional."  The  first 
step  in  calculating  the  channel  current  is  to  characterize  the  relationship  between  the 
applied  biases  and  the  surface  potential  and  inversion  charge,  including  short- 
channel  effects  based  on  two-dimensional  analysis  through  the  use  of  Poisson's 
equation.  Characterization  of  the  effective  quasi-two-dimensional  channel  length, 
based  on  two-dimensional  analyses  near  the  source  and  drain  junctions  for  arbitrary 
bias,  is  part  of  the  model  algorithm.  The  corresponding  terminal  quasi-static  charge 
modeling,  yielding  C^  continuity,  is  then  developed  using  the  DC  channel-current 

solution. 

3.2  Short-Channel  Effects 

In  scaled  devices,  the  effects  of  source/drain  charge  sharing  and  drain- 
induced  effects  (i.e.,  DIBL  and  DICE)  must  be  accounted  for  in  the  model  formalism. 
More  than  twenty  years  ago,  Yau  [Yau74]  introduced  a  simple  model  for  the 
threshold  voltage  of  short-channel  MOSFETs  based  on  a  geometrical  partition  of  the 
depletion  charge  among  the  gate  (QB(eff))'  source  and  drain.  As  the  channel  length 
is  shortened,  the  charge  shared  by  the  source/drain  becomes  more  significant,  with 
the  effective  gate-controlled  depletion  charge  and  hence  the  threshold  voltage  being 
reduced.  Physically,  the  threshold  voltage  dependence  on  channel  length  is  governed 
by  the  two-dimensional  Poisson's  equation  in  the  body. 


dx2  dy2  e. 
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which  characterizes  how  the  (lateral)  source/drain  junction  fields  affect  the  potential, 
inversion  charge,  and  even  the  field  in  the  oxide  [Vis85].  Accounting  for  charge 
sharing  in  compact  MOSFET  models  has  evolved  via  improved  characterization  of 
threshold  voltage  for  the  past  two  decades.  These  improvements  were  initially  based 
on  better  geometrical  charge  partitions  [Tay78],  and  more  recently  on  quasi-two- 
dimensional  analyses  [Kha93,  Liu93].  However,  all  of  the  models  are  limited  by  the 
assumption  that  charge  sharing  is  not  a  function  of  the  gate  bias,  and  their 
applicability  is  limited  to  characterizing  threshold  voltage  only. 

Many  models  are  based  on  the  gradual-channel  approximation  (GCA), 

d^\i/     d  \i/ 
which  assumes  that  — ^» — t  =  0,  implying  that  there  are  no  significant  two- 
dx''     dy^ 

dimensional  effects  in  the  channel  region  of  the  device.  For  characterizing  the  short- 
channel  effects,  a  new  assumption. 


^>^.C.f(x,y),  (3.2) 

dx''      dy^ 


is  proposed  to  be  valid  over  some  effective,  or  quasi-two-dimensional  region  of 
length  Lq2d;  C  is  spatially  constant  but  bias-dependent.  As  shown  in  Figure  3.1  (for 
an  n-channel  device),  the  lateral  field  associated  with  the  source  and  drain  junctions 
produces  regions,  with  lengths  L^  and  Lj  respectively,  in  which  (3.2)  can  not  be 
valid.  In  order  to  use  (3.2),  the  quasi-two-dimensional  channel  length  must  be 
defined,  using  two-dimensional  analyses  to  characterize  L^  and  Lj: 

LQ2D  =  L-L,-Ld  (3.3) 
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where  L  is  the  metallurgical  channel  length.  Such  analyses  have  been  done 
previously  [Gre93b],  [E1M77],  but  not  with  the  generality  concerning  bias  and 
scalability  needed  here.  The  definition  of  Lq2d  and  its  bias  dependence,  which  are 
given  later,  validates  (3.2)  and  allows  simpler  accounting  for  the  two-dimensional 
effects  in  the  quasi-two-dimensional  channel. 

3.2.1  Quasi-Tvro-Dimensional  Potential 

Finding  a  solution  to  the  two-dimensional  Poisson's  equation  in  (3.1)  is 
accomplished  by  assuming  that  the  two-dimensional  potential  can  be  written  as  the 
linear  combination  of  the  one-dimensional  solution  and  a  two-dimensional 
perturbation: 

VCx.y)  =  Vo(x) -H  8v|;(x,  y)  (3.4) 

where  \|/„  is  the  one-dimensional  solution  which  is  not  a  function  of  y  and  By 

represents  the  two-dimensional  perturbation  due  to  the  effects  of  charge-sharing  and/ 
or  drain  bias.  With  (3.4),  (3.1)  is  transformed  to 


d(Vo  +  5\|/)     d^gm         Po  +  5p 

T^ +  ^-T^  =  — •  (3-5) 

dx^  dy2  £ 


For  weak  inversion  where  the  inversion  charge  is  negligible,  8p  =  0 ,  and  since  the 
one-dimensional  Poisson's  equation  is 


^-    ^°  r3fi^ 
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(3.5)  yields  a  two-dimensional  Laplace's  equation  for  the  perturbation  in  the 
potential: 

#  +  #  =  0-  (3.7) 

For  low  (but  not  negligible)  longitudinal  fields,  (3.7)  can  be  solved  by 
quasi-separation  of  variables: 


dx-'  dy^ 


which  is  consistent  with  (3.2);  ^  models  the  quasi-two-dimensional  effects  on  \|/. 
Integrating  (3.8)  twice  in  y  over  the  quasi-two-dimensional  channel  length  results  in 
a  parabolic  potential  perturbation  8\|/(y),  which  yields  an  expression  for  C,: 


^_^Svt/(LQ,p)-8v(0)^^5Ey(0)  ^^^^ 

Lq2D  ^Q2D 


where  the  first  term  predominantly  represents  the  drain-induced  effects  and  the 
second  term  predominantly  represents  the  charge-sharing  effects.  Under  "high"  drain 
bias  where  the  drain-induced  effects  tend  to  dominate  the  two-dimensional  effects, 
the  parabolic  potential  has  a  minimum  near  the  source  junction,  which  implies  the 
approximations   [5\);(Lq2p)  -  5^/(0)] ->  Vps   and  Lq2I3 ->  L   in  (3.9).  Underlying 

these     approximations     are     the     boundary     values,     6\|/(-L5)  =  <t>s-Vo     and 

8\|/(Lq2[) -H  Lj)  =  <|)s  ~  Vo  +  ^Ds '    where    <^^    is    the    source    potential,    i.e.,    the 
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equilibrium  potential  barrier  of  the  source/drain  junction,  ([ij  =  -5  +  V,ln[  —  j  where 

^q  \  nj  y 

Ng  is  the  channel  doping  density.  The  potential  \\l  is  referenced  to  the  unbiased 
neutral  low-doped  body.  Using  these  two  approximations  results  in 

C  =  2-f  +  2--^.  (3.10) 

The  second  term  in  (3.10),  which  dominates  at  "low"  drain  bias,  can  not 
be  explicitly  defined.  In  order  to  estimate  it,  let  V^j  =  0  and  integrate  (3.8)  twice 
in  y  over  the  metallurgical  channel  length.  The  resulting  parabolic  5l|r(y),  with  the 


54)=0, 


assumption  5\|;  -    =  0 ,  yields 


5V-L,.2^--I^.  (3.11) 


The    parabolic    potential    assumption    defines    a   linear    variation    for    the    field 
perturbation.  Hence 

5Ey(0)s5Ey(-L,)-CL,(l+X),  (3.12) 

where   X  >  0    accounts  for  the  uncertainty   introduced  by  assuming  a  parabolic 
potential  in  the  strongly  two-dimensional  source  region.  For  Vps  =  0,  ^s2 — ^ — 

Lq2D 

in  (3.12),  which  gives  : 
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5E  (0)  8E  (-L,) 

-^ —  -  ^-  (3.13) 


Lq2D  Lq2]3  +  2L5(  1  +  X,) 

Due  to  the  symmetry,  L,.  =  Lj  and  Lqjd  +  2Lj  =  L,  which  transforms  (3.13)  to 


8Ey(0)  ^  SEyC-L^) 

^Q2D  ^  ■*"  ^LjA, 


Using  (3.14)  to  model  the  charge-sharing  effects  would  require  an 
iterative  solution  since  Lj  is  a  strong  function  of  ^ .  In  order  to  make  the  analysis  for 
the  two-dimensional  effects  suitable  for  compact  modeling,  the  approximation, 
L  +  2h^X  s  L(  1  +  FCS) ,  where  PCS  is  a  model  parameter  greater  than,  but  near  0,  is 
introduced  to  allow  a  closed-form  solution.  Then  (3.14)  with  (3.11)  becomes 

8Ey(0)  (|)s-Vo(x) 

Lq2d  L2(1+FCS)-  ^^•'^' 

However,  with  (3.15)  the  expression  for  t,  in  (3.10)  is  a  function  of  x  which  is 
inconsistent  with  the  model  in  (3.2).  Therefore  \|/„(x)  in  (3.15)  is  replaced  by  an 
average  value  of  the  (parabolic)  potential  across  the  one-dimensional  depletion 
region  (of  depth  \^„),  which  will  be  defined  later.  In  order  to  define  this  value,  the 
potential  is  integrated  over  the  depletion  depth  using  the  depletion  approximation, 
and  the  mean  value  theorem  is  applied:  '     * 
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Vo  =  —  j  'VoM'^^ 


Vs„-Vb 


(3.16) 


Using  (3.16)  in  (3.15)  now  yields  an  expression  for  C,  in  (3.10)  which  is  spatially 
constant: 


?^2. 


'ds  ,  o      ^s  -  Vo 


L2         L2(1+FCS)' 


(3.17) 


In  order  to  extend  the  model  for  strong  and  moderate  inversion,  it  is 
assumed  that  the  conditions  in  the  channel  at  threshold  adequately  represent  the 
conditions  below  the  inversion  layer  in  strong  inversion  due  to  the  screening  of  the 
depletion  charge  from  the  gate  by  the  inversion  layer.  This  is  effected  by  limiting  the 


one-dimensional  surface  potential  in  (3.16)  to  less  than  <|)^  =  2(|)p  =  V,ln 


(K^ 


which 


represents  the  onset  of  inversion.  Analogous  to  a  function  presented  by  McAndrew 
etal.  [MeA91],  ,   ; 


V,; 


In 


1  +e 


^1 


(3.18) 
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where  A  is  a  hard-wired  parameter  that  controls  the  abruptness  of  the  transition, 
replaces  \|/j^  in  (3.16).  The  use  of  (3.18)  smooths  the  transition  between  the  physical 
approximations  in  weak  and  strong  inversion;  such  intrinsic  smoothing  is  physical, 
in  contrast  to  single-function  models  [Aro94]  that  extrinsically  and  empirically  yield 
smooth  terminal  characteristics. 

3.2.2  Quasi-Two-Dimensional  Channel  Charge  Perturbation 

Based  on  the  characterization  of  C,,  the  inversion  charge  in  the  quasi-two- 
dimensional  channel  region  is 

Qc  =  Qco  +  5Qc  (319) 

for  all  bias  conditions.  The  perturbation  5Q;. ,  which  corresponds  to  8v|/,  is  expressed 
by  first  integrating  (3.8)  in  x  over  the  one-dimensional  depletion  region  (xj^)  to 
describe  the  perturbation  in  the  vertical  field  below  the  inversion  layer: 


"do  2 


Gauss's  law  is  then  applied  over  0<  x<Xj,  where  x  =  0  lies  just  above  the  Si-SiOj 
interface  and  Xj  is  the  inversion  layer  thickness: 

6Q,  =  e35E,-e,8E,.  (3.21) 

Assuming  vertical  fields  in  the  oxide  and  substituting  (3.20)  into  (3.21)  yields  an 
expression  for  quasi-two-dimensional  channel  charge  perturbation: 
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8Qc  =  (c„;  +  ;^)6v,-'-i^°,  (3.22) 

where  C„j'  represents  the  variation  in  the  polysilicon  potential  along  the  channel  as 
defined  later. 

The  two-dimensional  effects  (Q  have  now  been  approximated  for  all  bias 
conditions,  and  the  corresponding  perturbations  of  surface  potential  and  inversion 
charge  were  described.  However,  the  total  solution  depends  on  the  one-dimensional 
solution  (Vo' Vso' Qco)  •  ^he  solution  to  the  one-dimensional  Poisson's  equation  is 
now  presented. 

3.3  One-Dimensional  Surface  Potential  and  Inversion  Cliarge 

In  this  section,  the  one-dimensional  relationships  among  the  surface 
potential  i^^„)  the  inversion  charge  (Q^.^)  and  the  intrinsic  biases  (Vqj,  Vgs  and 
Vds)  in  the  quasi-two-dimensional  channel  region  are  derived  in  order  to 
characterize  the  channel  current. 

Near  the  center  of  the  n-channel  MOSFET  structure  in  Figure  3.1, 
Poisson's  equation  is  applied  in  one  dimension  accounting  for  hole,  depletion  and 
electron  charge: 


d^  =  _5!!B^[e-P(¥.-VBs)_i_eP('Co-W  +  eP(VBs-*.)]  (323) 

dx^  e. 
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1  kT 

where    Ng(x)     is    the    channel    doping    concentration,     P  -  77 .     V,  =  —     and 


(!.b  =  2(t)psV,ln 


fNgl  dv„ 

Multiplying  both  sides  of  (3.23)  by  2— —  and  integrating  over 


\|/j,  results  in  a  relationship  between  the  electric  field  and  potential  at  the  surface 
(x  =  0): 

Bl  =  QJe-P«^.»-VB='-l+P(X|/,„-VBs)  +  eP<v„-'t'.)_eP(VBs-0.)]  (3.24) 


2qNBV,  _ 

where    Q^  = and    Ng,    discussed    later,    represent    an    effective   doping 


concentration  based  on  the  vertical  doping  transform  [Aro87]  of  the  retrograde 
profile;  the  vertical  field  deep  in  the  bulk  is  implicitly  assumed  to  be  zero. 

The  term  on  the  left-hand  side  of  (3.24)  is  the  surface  field  in  the 
semiconductor,  which  with  the  assumption  that  the  field  in  the  oxide  is  one- 
dimensional,  can  be  expressed  via  Gauss's  law  in  terms  of  \\l^g  and  \|/  ,  the  potential 
drop  across  the  polysilicon-gate,  which  will  be  analyzed  later,  as 


Eso  =  -f  (Vgs  -  VpB  -  Vso  -  V)  (3-25) 


^ox  Qf 

where  C^,^  =  —  is  the  gate  oxide  capacitance,  Vpg  =  <t>^3  -  -;—   is  the  flatband 


voltage,  0^3  =  -  ;r^  -  V,ln| —    is  the  gate-body  work-function  difference  and  Qf 


El 

2q 
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represents  the  fixed  charge  (areal)  density  in  the  gate  oxide.  With  (3.25),  (3.24)  can 
be  solved  for  \^f^g  in  a  few  iterations  using  a  second-order  Newton's  method.  Note 
that  (3.24)  is  identical  to  the  classical  charge-sheet  equation  [Bre78]. 

3.3.1  Nonuniform  Vertical  Doping  Profiles 

The  vertical-doping  transform  [Aro87]  enables  a  step  junction  to  be 
represented  with  an  effective  constant  doping.  Although  the  analysis  was  originally 
developed  to  approximate  a  Gaussian  profile  with  an  equivalent  (bias-dependent) 
uniform  doping  concentration,  the  analysis  is  also  valid  for  the  case  of  retrograde 
channel  doping  profiles.  This  transform,  for  an  arbitrary  vertical  doping  profile,  is 
based  on  four  conditions: 

•  ''do>''LH 

•  Conservation  of  depletion  charge  under  the  gate 

•  Constant  surface  potential  \\i^g 

•  Continuous  electric  field  E^  at  x  =  x^h 
which  lead  to  the  bias-dependent  effective  doping  given  by 


N„  =  N„ 


Vi 


B    -   '"B,,,      _v       1  *    ■     N 


N„ 


1  + 

BSl  '"B 


1  + 


NefYs. 
NrhI 


Vi 


1-1 


(3.26) 


where  Ng  is  the  surface  doping  concentration,  N^^^  is  the  doping  concentration  in 


the  p*  layer  and  \j/,  = 


2e. 


is  the  built-in  potential  drop  across  the  lightly  doped 


layer. 
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However,  this  analysis  assumes  that  the  lightly  doped  layer  is  always 
depleted,  which  is  not  the  case  as  bias  conditions  approach  flatband.  The  term 

Vso-Vbs 


Vi 


is  therefore  replaced  in  (3.26)  with  a  smoothly  varying  function 


In 
r  =  1  +  — 


(3.27) 


Vso~^BS 

which  is  approximately  when  the  lightly  doped  region  is  depleted  and 

approximately  1  when  it  is  not  fully  depleted.  Thus,  Ng  — >  Ng  as  the  bias  conditions 

approach  flatband,  which  is  shown  in  Figure  3.2  for  a  typical  device.  The  smoothing 
function  guarantees  the  continuity  of  the  effective  doping  concentration  and  hence 

¥,0- 

The  analysis  of  the  effective  doping  also  results  in  an  expression  for  the 
depletion  depth.  Starting  from  the  assumption  that  the  surface  region  is  depleted,  the 
(maximum)  depletion  depth  is 


'-"  =  XLHJl+i;T^(r-I).  (3.28) 


However,  when  the  surface  region  is  not  fully  depleted  the  depletion  depth  is  given 
by  the  depletion  approximation:  •'i   :  .  .  . 
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Figure  3.2  The  Vertical  Doping  Transform  for  a  typical  device. 
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qNfi 


(3.29) 


Another  smoothing  function  is  then  used  to  find  the  depletion  depth  under  all  bias 
conditions: 


B   I 

1+e  '^     W 


ln[l+eB] 


(3.30) 


3.3.2  Nonuniform  Lateral  Doping  Profiles 

In  scaled  devices,  it  is  also  necessary  to  account  for  nonuniform  lateral 
doping  profiles  due  to  halo  implants  [Tau97]  and  also  the  Reverse  Short-Channel 
Effect  (RSCE)  [OrI87].  A  preliminary  accounting  of  these  effects  is  done  by  using  a 
lateral-doping  transform  which  assumes  quasi-Gaussian  profiles  (with  peak  doping 
of  No  and  a  characteristic  length  y,,)  near  the  source/drain  junctions  to  modify  the 
surface  doping  concentration  in  (3.26)  as  a  function  of  the  channel  length  [Bru95]: 


Nb  =  Nb,l- 


+  2^° 


,-J)" 


(3.31) 


Thus,  the  effective  Ng  is  higher  than  the  background  doping  density,  Ng.^^^),  as 
illustrated  in  Figure  3.3 
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Figure  3.3  The  Lateral-Doping  Transform  for  a  typical  device. 
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3.3.3  Polysilicon-Gate  Depletion 

In  scaled  technologies,  the  use  of  dual  polysilicon  gates  is  predominant 
due  to  the  difficulty  of  designing  deep-submicron  buried  channel  PMOS  devices 
[Sha96]  and  the  simplicity  of  simultaneously  implanting  the  polysilicon-gate  and 
source/drain  junctions  [Hua93].  However,  the  implant  processing  can  limit  the 
electrically  active  doping  concentration  in  the  polysilicon  [Sch93],  which  can  result 
in  a  modification  in  the  device  characteristics  (both  DC  and  AC)  due  to  the  potential 
drop  across  the  depletion  region  formed  in  the  polysilicon  gate.  Although  it  is 
possible  to  achieve  good  fits  to  DC  data  without  modeling  the  polysilicon-gate 
depletion  effect  [Aro95],  significant  errors  in  the  capacitance  characteristics  will 
lead  to  erroneous  transient  and  AC  circuit  simulations.  It  is  therefore  necessary  to 
account  for  this  effect  in  order  to  model  scaled  devices. 

Assuming  that  the  two-dimensional  effects  in  the  gate  oxide  and 
polysilicon  are  negligible  and  that  the  doping  profile  is  uniform  near  the  poly-Si02 
interface  [Hua93],  Gauss's  law  and  the  depletion  approximation  give  the  potential 
drop  across  the  gate  material,  based  on  the  electric  field  in  the  polysilicon  at  the 
interface: 


Epo  =  -f^CVos  -  <1'ms  -  Vso  -  ¥po)  (3.32) 


and 


Vpo  =  ;T7^^^E2  ■,     rv...,^    .  (3.33) 
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where  NpQ^Y  is  the  active  doping  concentration  in  the  polysilicon  gate.  Substituting 
(3.33)  into  (3.32)  results  in  a  quadratic  equation  for  E     , 


2^^f^Ep\  +  ^Ep„  -  (Vgs  -  «I>Ms  -  Vso)  =  0 ,  (3.34) 


which  can  be  solved  analytically  and  substituted  into  (3.33)  to  find  the  potential  drop 
across  the  polysilicon  gate. 

3.3.4  Quantum  Mechanical  Effects  in  the  Inversion  and  Accumulation  Layers 

In  extremely  scaled  technologies,  the  high  fields  at  the  silicon  surface 
produced  by  thinner  gate  oxides  and  higher  channel  doping  lead  to  a  splitting  of  the 
conduction  and  valence  bands  into  discrete  energy  sub-bands  through  quantum 
mechanical  effects.  Although  this  effect  can  be  modeled  by  assuming  an  "electrical" 
tgx  larger  than  the  physical  gate  oxide  thickness  [Aro95],  this  introduces  the  problem 
of  non-physical  model  parameters  and  poor  capacitance  characteristics  in  the 
accumulation  and  moderate  inversion  regions  of  operation.  The  quantization  of  the 
carriers  in  the  inversion  and  accumulation  layers  can  be  treated  as  an  effective 
bandgap  widening  [Dor94]  through  a  correction  to  the  intrinsic  carrier  concentration: 

npM  =  n,exp[-^].  ■•  (3.35) 

With  the  assumption  of  a  triangular  potential  well  and  the  use  of  a  variational  method 
taking  into  account  only  the  first  sub-band,  the  effective  change  in  the  bandgap  is 
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^^s  =  K4rf^-'   •  ^^-^^^ 


where  K  =  4.1x10"  eV-m  was  extracted  from  threshold  voltage  measurements  on 
devices  with  uniformly  doped  substrates  and  14nm  gate  oxides.  In  order  to  account 
for  increased  quantum  mechanical  effects  in  thinner  oxides  where  k  may  be  larger, 
and  also  to  turn  off  the  quantum  mechanical  effects  in  UFET,  (3.36)  is  modified: 


AE^  =  QM^k(A,e2„)'"  (3.37) 


where  QM  is  a  model  parameter. 

Accounting  for  quantum  mechanical  effects  in  the  UFET  model  formalism 
is  done  through  a  modification  of  Poisson's  equation.  As  proposed  by  Arora  et  al. 
[Aro95],  the  spatial  dependence  of  the  quantum  mechanical  effect  is  ignored,  and  the 
resulting  modification  of  (3.24), 


Ei  =  Qr[^e  _i  +  p(^^^_V3s)  +  -^e  -eP(^=s-».)J ,  (3.38) 


can  still  be  solved  in  a  few  iterations  although  the  nonlinearity  of  the  equation  has 

npM 

increased.  Note  in  (3.38), >  1  and  (3.38)  reverts  to  the  classical  expression  for 

n. 


AE„^0. 
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3.3.5  Inversion  Charge 

After  solving  (3.38),  the  magnitude  of  the  one-dimensional  inversion 
charge  can  be  obtained  by  first  integrating  Poisson's  equation,  following  the 
derivation  of  (3.24),  over  the  predominate  inversion  layer  to  define  the  vertical  field 
below  the  inversion  charge,  E^^; 


E2„-E2„  =  QJ-^e  -eP^^Bs-P.'!.  (3.39) 


Then,  Gauss's  law  is  used  to  define  the  inversion  charge: 


Qco(y)  =  e. 


pQM    PCv,,,- 


'b) 


'E.^„-Q{^e  -eP(VB.-wJ_E^J.  (3,40) 


3.4  Quasi-Two-Dimensional  Inversion  Charge 

In  order  to  find  the  magnitude  of  the  two-dimensional  inversion  charge  in 
the  quasi-two-dimensional  channel,  Poisson's  equation  is  used  in  conjunction  with 
Gauss's  law,  as  was  done  in  the  derivation  of  (3.40): 


Qc(y)  =  e, 


L  n,  J        ^J  (3.4!) 

where  Vy  is  the  additional  quasi-Fermi-level  separation  at  y  implied  by  the  drain 
bias.  Although  (3.41)  is  of  the  same  form  as  (3.40),  this  is  consistent  with  the  charge- 
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sheet  approximation  where  the  inversion  layer  has  an  infinitesimal  thickness  and  the 
contribution  to  Gauss's  law  due  to  the  two-dimensional  effects  vanishes: 

e,Cxi-^0  (3.42) 

where  x,  is  the  inversion  layer  thickness. 

Note  that  in  (3.41)  E^  is  no  longer  the  one-dimensional  surface  field,  but 

may  be  influenced  by  the  two-dimensional  effects.  The  two-dimensional  surface 
field  in  the  semiconductor,  with  the  assumption  that  the  field  in  the  oxide  is  one- 
dimensional,  can  be  expressed  via  Gauss's  law  in  terms  of  \\i^g  +  b\\i^  as 


Es  =  Y^{  Vgs  -  VpB  -  Vso  -  8Vs  -  Vpo) .  (3-43) 


where  the  use  of  \|/     assumes  that  the  polysilicon-depletion  effects  are  not  perturbed 

by  the  two-dimensional  effects.  It  is  also  assumed  that  (3.35)  and  (3.37)  adequately 
represent  the  quantum-mechanical  effects  using  the  one-dimensional  surface  field. 

Using  (3.40)  and  (3.41)  it  is  now  possible  to  define  the  perturbation  in  the 
two-dimensional  channel  charge  as 

6Q,  =  Q,-Q„.  (3.44) 

This  result  can  then  be  used  with  (3.22)  to  fully  characterize  the  two-dimensional 
effects  on  potential  and  inversion  charge.  The  combination  of  (3.22) 
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(3.44)  must  be  solved  iteratively  for  5\|/j,  (~4  iterations  worst-case)  using  a  Newton- 
Raphson  algorithm.  Note  that  in  weak  inversion  the  two-dimensional  effects  will 

result  in  a  5v(/jS while  in  strong  inversion  the  resulting  effect  is  a 

5Qj.s — ^      °.  The  iterative  solution,  in  terms  of  the  intrinsic  biases,  guarantees  a 

smooth  physical  transition  between  these  two  effects  across  the  moderate-inversion 
region. 

3.5  Channel  Current  Model 

The  current  at  any  position  along  the  channel  can  be  expressed  as  the  sum 
of  the  drift  and  diffusion  components: 


I  =  WnV,-^-WQ,v  (3.45) 


where  |i  is  the  carrier  mobility  and  v  is  the  electron  drift  velocity. Using  our  proposed 
spatially  constant  two-dimensional  effects  assumption  in  (3.45)  enables  analytic 
accounting  for  all  of  the  two-dimensional  effects  at  any  point  along  the  quasi-two- 
dimensional  channel,  and  hence  a  solution  for  the  drift/diffusion  channel  current.  In 
order  to  find  the  channel  current,  (3.45)  must  be  integrated  along  the  channel  to  find 
a  closed-form  solution,  which  is  needed  in  a  compact  model  to  minimize  computation 
time.  In  order  to  achieve  this  goal,  the  choices  for  the  mobility-(vertical)  field  model, 
the  velocity  saturation  model  and  inversion  charge  variation  model  in  UFET  have 
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been  influenced.  The  need  for  an  explicit  channel  current  equation  necessitates  the 
forms  of  the  models  to  be  relatively  simple  and  easily  integrated  along  the  channel. 

3.5.1  Mobility-Vertical  Field  Model 

The  transverse  field  effect  on  mobility  is  dependent  on  E^  in  the  channel; 

it  is  modeled  in  terms  of  the  average  Ejj(y)  along  the  channel  [Ong87]: 


Heff  =  =^ =5—.  (3.46) 

l+0,E,(y)  +  G2E,^(y) 


where  n„  is  the  bulk  mobility  and  0|  and  02  are  empirical  constants.  The  average 
transverse  field  is  defined  as 

E"x(y)  =  -^[VGs-VFB-¥s(y)-Vp(y)]  +  ^.  0.47) 

3.5.2  Velocity-Lateral  Field  Model 

For  large  drain  bias  conditions,  the  high  longitudinal  electric  fields  in  the 
channel  produce  a  nonlinear  carrier  velocity-field  relationship.  At  sufficiently  high 

longitudinal  fields,  the  carrier  velocity  saturates  (v^j,,  -  8x10     cm/s  for  electrons). 
The  carrier  velocity  is  modeled  by  [Aro94] 


^=ffd7 
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where  ji^fj  is  the  low-(longitudinal)-field  mobility  and  T  is  defined  based  on  an 
average  (longitudinal)  electric  field  as 

T  =       ^^^'_   .  (3.49) 


1.5  + 


\lEy 


3.5.3  Channel  Current 

With  the  mobility  implied  by  (3.48),  the  drift/diffusion  current  in  (3.45)  is 
expressed  as  ■      ' 

^^ff"d7 

l+T 2L 


Note  that  although  the  mobility  implied  by  (3.48)  is  only  strictly  valid  for  the  drift 
current  term,  it  is  used  for  the  diffusion  term  also.  However,  since  the  longitudinal 
electric  field  under  bias  conditions  of  predominant  diffusion  current  is  small,  the 
error  introduced  is  negligible. 

Although  the  inversion  charge  along  the  quasi-two-dimensional  channel 
can  be  defined  in  terms  of  the  surface  potential  and  quasi-Fermi  level  separation  with 
(3.41),  the  variation  of  the  quasi-Fermi  level  separation  along  the  channel  is  not 
known,  and  the  inversion  charge  must  therefore  be  expressed  solely  in  terms  of  the 
surface  potential  in  order  to  find  a  closed-form  solution  for  the  channel  current.  For 
convenience,  a  reference  point,  i.e.,  the  virtual  source,  is  defined  in  the  channel  (at 
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y  =  0),  and  the  potential  and  inversion  charge  at  all  points  along  the  channel  are 
referenced  to  this  point.  With  the  solution  for  the  two-dimensional  potential  and 
inversion  charge  at  y  =  0,  Vs(0)  and  Q(,(0)  respectively,  the  change  in  inversion 
charge  along  the  channel  is  defined  via  Gauss's  law  as 

AQc(y)^-esAEs-es^Qb.  (3.51) 

where  AE,(y)  ^E,(y)-E,(0)  and  AQ^(y)  ^Q^Cy) -Q,(0). 
The  change  in  the  surface  field  can  be  modeled  as 

Cox 
AEs(y)  =  -^[Av,(y)  +  Av|/p(y)] .  (3.52) 

s 

The  earlier  analysis  for  the  one-dimensional  polysilicon  depletion  potential  results  in 
an  expression  for  the  polysilicon  depletion  potential  along  the  channel, 


Vp(y)  =  Vp(0)  +  A\|/p(y) 

=  V  +  ^Vp(y) 


(3.53) 


[E     +AE,(y)]2 


which  can  be  solved  to  find  Ai);  (y) : 


_  eJ2Ep„AE,(y)-HAE,(y)] 


AVp(y)  =     -      ""      -' ^— .  (3.54) 


Substituting  (3.54)  into  (3.52)  results  in 
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Av.(y).'-^''^-'^^'^'"^'^^^^^' 


2qNp 


(3.55) 


In  the  resulting  quadratic  equation  in  AE^  only  the  first-order  term  is  kept 

in  order  to  allow  the  resulting  channel  current  equation  to  be  evaluated  analytically 
in  the  calculation  of  the  drain  saturation  potential.  This  is  equivalent  to  assuming  a 
constant  depletion  capacitance  along  the  channel  in  the  polysilicon  gate.  The 
variation  in  the  field  is  now  expressed  as 


Y^Av,(y) 
AE,(y)  =  -- ^-p (3.56) 

1  + X„ 


£  E 

where  x    =     ^   ''°    accounts  for  the  depletion  in  the  polysilicon-gate. 

The  variation  in  the  depletion  charge  along  the  channel  is  modeled  using 
the  depletion  approximation: 


AQb(y)  =  J2qe^N^^^,^^-J2qe^Ng[^^,^„  +  A^v,iy]].  (3.57) 

However  in  order  to  obtain  the  desired  analytical  solution  of  the  drain  saturation 
potential,  a  Taylor  expansion  of  (3.57)  is  taken  about  Av|i5(y)  =  0  and  the  first-order 
term  is  kept: 
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AQb(y)  =  -^Av,(y). 


(3.58) 


This  is  equivalent  to  using  a  constant  depletion  capacitance  to  model  the  variation  in 
depletion  charge.  The  variation  in  the  inversion  charge  density  along  the  channel  is 
now  expressed  as 


AQ,(y)  = 


AVs(y) 


(3.59) 


=  (Cox'  +  Cd)Av|/,(y) 


where  C^^'  in  (3.59)  accounts  for  the  variation  in  the  polysilicon  potential  along  the 
channel. 

From  the  integration  of  (3.50)  over  the  length  (Lq2d)  of  the  quasi-two- 
dimensional  channel,  the  steady-state  channel  current  is  now  modeled  as 


I      -    ^  n~[viOfL      ^    nrmi  +  5?(^hQ^liQ2D) 

^Q2D        [  ^  2(C„^  +  Cj) 


(3.60) 


where 


^l.ff(0) 


^cSf  = 


l-BiAy^(LQ2D)-B2Ay2(LQ2p) 


^leff(0) 


1+T 


1-B|Av|/,(Lq2p)-B;Av|/2(Lq;d) 


AV|/,(Lg2o) 


(3.61) 


\4^< 


Heff(O)  =  = °- =1— '  (3.62) 

l+e,E,(O)  +  02E,2(O) 


[0,+ 2626,(0)]- 


B,  =  = %~ (3.63) 

(fB  +  i)[i+e,E,(0)  +  e2E,2(0)] 


and 


<-^''' 


(3.64) 


(2fB  +  l)[l+0,E,(O)  +  e2E/(O)] 

Note       that       to       obtain       (3.63)       and       (3.64),       the       approximation 

^Vs(y)  =  ^Vs(LQ2D)f[7— J"  was  used  to  evaluate    J  Ay^dy  and    jAv52dy;this 
'^^^  0  0 

approximation  is  not  used  anywhere  else  in  the  model,  and  in  no  way  implies  that  a 
fixed  dependence  of  Av|/5(y)  governs  (3.60). 

Together  (3.40)  and  (3.59)  give  the  inversion  charges  at  both  ends  of  the 
channel,  which  are  used  in  (3.60)  to  define  the  channel  current.  The  only  unknowns 
in  (3.60)  are  Lq2d  and  A\|/5(Lq2p).  The  next  step  is  then  to  define  the  physical 
locations  of  the  y  =  0  and  y  =  Lqjp  points  along  the  quasi-two-dimensional 
channel,  which  leads  to  an  expression  for  A\|;5(Lq2c)  and  the  primary  Vps 
dependence. 
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3.6  Two-Dimensional  Source-Side  Analysis 

Focusing  on  the  source  region  of  the  MOSFET,  Gauss's  law  is  applied  to 
an  incremental  portion  (dy)  of  the  region  defined  by  -L^  <  y  <  0  and  0<x<Xjj,, 

where  x  =  0  lies  just  above  the  Si  -  SiOj  interface  and  Xj^  is  the  depletion  depth  at 

the  virtual  source,  which  is  used  as  the  fixed  lower  boundary  of  the  Gaussian  surface. 
Application  of  Gauss's  law  results  in 

Qc(y)  +  Qb(0)  +  Qf  =  esEx(''do.  y)  -  eoxE„x(0,  y )  -  e J  ^dx .  (3.65) 


,2 

With  the  assumption  of  — ^^^^fCx),  the  vertical  field  at  x  =  x^^  can  be 
found  by  integrating  the  two-dimensional  Poisson's  equation, 


dV  _  INb     dVs 


twice  over  the  one-dimensional  depletion  depth  and  solving  for  E^{x^^,  y) : 

Ex(Xdo'y)  =  — - 

''do 


l^.       dy^j2-  (^■^') 


Now  using  the  solution  at  the  virtual  source, 

Qc(0)  +  Qb(0)  +  Qf  =  e,E,(x,„,  0)  -  e„E„,(0,  0)  -  e,x,„? ,  (3,68) 
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which  provides  a  convenient  reference  point,  (3.65)  can  be  written,  after  subtracting 
(3.68),  as 

^^  =  ^i^^l,^-,,0)J-.S.M-».,, 

which  is  similar  to  the  analysis  of  current  saturation  in  strong  inversion  [E1M77].  In 
that  work,  the  channel  was  segmented  into  a  region  near  the  source  where  the 
gradual-channel  approximation  can  be  applied  and  a  high-field  region  near  the  drain 
where  the  carriers  have  a  saturated  drift  velocity.  Since  the  current  is  constant  along 
the  channel,  the  saturated  drift  velocity  implies  that  the  carrier  concentration  is  a 
constant  in  this  region,  and  there  is  therefore  no  term  accounting  for  any  variation  in 
inversion  charge  Qj(y),  in  contrast  to  (3.69). 

However,  it  should  be  noted  that  C,  given  by  (3.17)  is  much  smaller  than 

.  .      ■     ■       ■  J''f", :_ 

- — —  for  long-channel  devices.  The  ^  implied  by  (3.2)  is 
dx^ 


d  w.„ 
C  =  f-^  (3.70) 


where  f  «  1  .  This  value  is  used  as  a  minimum  value  for  ^  in  the  channel-length 
modulation  analyses: 


d  VI, „ 
C^C  +  f--^.  (3.71) 

dx'' 


^k; 
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The  value  of  f  in  (3.71)  is  not  critical,  and  the  same  results  are  obtained  for  a  wide 
range  of  values  ( 10-2  <  f  <  10"^) . 

In  (3.69),  the  two-dimensional  Q;.(y),  characterized  by  (3.41),  is  not 
expressed  in  a  manner  which  allows  a  closed-form  solution  to  be  found.  To  simplify 
the  model  for  L^,  an  approximate  inversion  charge  distribution  which  is  valid  near 
the  source  junction  is  therefore  derived.  With  the  assumptions  that  the  drift  and 
diffusion  components  of  the  current  in  (3.45)  are  in  approximate  detailed  balance  in 
this  two-dimensional  region  and  that  the  low  lateral  electric  field  mobility  should  be 
used  in  the  two-dimensional  source  region  [Lun96],  (3.45)  then  gives 


'dy  -^'^dy 


(3.72) 


Integrating  (3.72)  from   y  =  0   to  an  arbitrary  position  then  yields  the  derived 
approximate  inversion  charge  distribution: 


>f.,-¥,(0) 


Q,(y)sQ,(0)e 


(3.73) 


Substituting  (3.73)  into  (3.69)  now  gives 


ay  Es^do  ^s'^do 


v.-v,(0) 


+  ^, 


(3.74) 


which  must  be  solved  to  characterize  L., 


Although  (3.74)  can  be  solved  numerically  via  Runge-Kutta  or  a  similar 
method,  the  exponential  term  accounting  for  the  inversion  charge  makes  it 
impossible  to  find  a  closed-form  solution.  Since  the  surface  potential  increases  on  the 
order  of  a  few  thermal  voltages  across  the  two-dimensional  region  in  strong  inversion 
where  the  characterization  of  the  inversion  charge  is  important,  a  first-order  Taylor 
expansion  of  the  exponential, 


.Qc(0) 


V.-WO)      n      _Q^(0)/v, 


b2^— -'[v^-vi/^CO)],  (3.75) 


can  be  used  to  enable  a  closed-form  solution  to  be  found  while  adequately  accounting 
for  the  inversion  charge.  Although  this  expansion  underestimates  the  inversion 
charge  in  the  two-dimensional  region  as  shown  in  Figure  3.4,  the  length  of  this  region 
is  controlled  mainly  by  the  potential  variation  that  occurs  near  the  quasi-two- 
dimensional  channel,  which  is  well  represented  by  (3.75).  The  figure  also  shows  that 
approximately  80%  of  the  length  of  the  two-dimensional  region  is  defined  by  an 
increase  in  surface  potential  of  half  a  thermal  voltage,  for  which  the  first-order 
Taylor  expansion  is  accurate. 

With  this  approximation,  the  differential  equation  (3.74)  can  now  be 
written  as 


-T-i  = — ^[Vs-Vs(0)]  +  C.  (3.76) 

"y  S^do 


which  can  be  solved  to  give  the  surface  potential  as  a  function  of  position: 
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0.5  1.0  1.5 

[Vs(y)-Vs(o)]/v. 


Figure  3.4  The  Taylor  expansion  of  the  inversion  charge  underestimates  the 
magnitude  of  the  charge. 
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Vs(y)  =  v)/.(0)  + 


dy 


^^^..sinhgj.^coshg)-!] 


(3.77) 


where  L 


2[Co;  +  Cj-Q^(0)/V,]- 
The  virtual  source  is  located  at  the  potential  minimum  along  the  channel, 

=  0  condition.  Substituting  this  condition  into 


which  is  equivalent  to  the  _!- 
dy 


y  =  o 


(3.77)  along  with  the  surface  potential  at  the  metallurgical  junction,  \|/j(-L5)  =  (t)^ , 
results  in  an  equation  which  relates  the  change  in  surface  potential  to  the  length  of 
the  two-dimensional  source  region: 


<)s  =  V3(0)  +  Cl2[cosh(^i^j-l] 


(3.78) 


This  equation  can  be  solved  analytically  using  the  identity  which  shows  that  the 


solution  to  A  =  Bsinh(x) +  Ccosh(x)  is  x  =  iJ^^^^-^^^^^^t^—^] 


B+C  j 


The  result  is 


U  =  Lin 


<t>s-Vs(0)^,_^  |[1>s-Vs(0)  ,  , 


?lc^ 


Clc^ 


(3.79) 


which  gives  the  length  of  the  two-dimensional  source  region  to  be  separated  out  from 
the  metallurgical  channel  length.  The  length  of  the  two-dimensional  source  region  is 
a  function  of  all  terminal  biases  (Vos,  Vgs  and  V^s)  as  well  as  channel  length  based 
on  the  dependencies  of  ^^{0),  Z,  and  QJO) . 


3.7  Two-Dimensionai  Drain-Side  Analysis 

Focusing  on  the  drain  region  of  the  MOSFET,  Gauss's  law  is  applied  to  an 
incremental  portion  (dy)  of  the  region  defined  by  Lgjp  ^  y  ^  Lqjo  +  L^  and 
0<x<Xjj(j,  where  x  =  0   lies  just  above  the   Si-SiOj   interface  and  Xjo  is  the 

depletion  depth  at  the  virtual  source  which  is  used  as  the  fixed  lower  boundary  of  the 
Gaussian  surface.  Analogous  to  the  previous  derivation  for  Lj,  application  of  Gauss's 

law  using  the  solution  at  y  =  Lqjo  as  the  reference  results  in 

dy^  es''do  ^s^do 

In  the  case  of  saturation,  the  current  analysis  in  the  quasi-two-dimensional 
channel  is  coupled  to  an  analysis  of  Lj  that  accounts  for  short-channel  effects,  as  well 
as  to  that  underlying  Lj.  In  the  two-dimensional  drain  region,  the  carrier  velocity  is 
saturated,  and  the  current  in  this  region  can  be  expressed  as 

IcH(sa,)  =  -WQc(y)V3a,.   '     '  (3-81) 

Since  both  the  (assumed  drift)  current  and  the  carrier  velocity  are  constant 
throughout  the  high-field  two-dimensional  drain  region,  the  inversion  charge  density 
must  also  be  constant;  this  charge  density  can  be  expressed  as 

Q,(y)  =  Q,(Lq2d).  (3.82) 
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With  the  constant  inversion  charge  density  assumption,  (3.80)  can  now  be 
expressed  as 

dV,       2(C    '  +  C.) 
dy^  Es^do 

The  derived  differential  equation  (3.83)  is  similar  to  the  classical  derivation  of 
velocity  saturated  channel  -length  modulation  [E1M77],  but  with  the  addition  of  the 
^  term  which  accounts  for  short-channel  effects. 

The  solution  to  (3.83),  with  the  boundary  conditions  for  the  potential 


A»j/  (Lq2I3  +  Lj)  =  Vp5    and  lateral  field   —1 

dy 


implied  by  (3.48),  yields 


y  =  L,,/^^^'"(l-Tf)tieff(LQ2D) 


Vds  -  Av|/s(Lq2d)  =  E,  JeSi"h(r )^  Cl?[cosh(^^)-  l]  (3.84) 


where  1,,  =    I     ^'^''°  Equations  (3.3),  (3.60),  (3.81)  and  (3.84)  provide  an 

implicit  solution  for  Lq2d  and  A\|/j,(Lq2[3)  ,  and  hence  IcH(sat)  • 

Although  the  channel  current  can  be  found  iteratively,  a  discontinuity  in 
the  derivatives  at  the  onset  of  velocity  saturation  would  result  due  to  the  assumption 
of  a  non-physically  abrupt  onset  of  velocity  saturation.  A  more  physical  solution, 
which  leads  to  a  C^,, -continuous  model,  is  found  using  a  system  of  equations  that  is 

asymptotically  correct  but  has  a  smooth  transition  between  the  linear  and  saturation 
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regions  of  operation.  The  first  step  in  finding  the  approximate  solution  is  to  assume 
that 


A¥s(Lq2d)  =  ^^!' 


1  +e 


ln[l  +eB] 


(3.85) 


B  is  once  again  a  hard-wired  constant  that  controls  the  abruptness  of  the  transition 
from  the  linear  to  the  saturation  region.  The  perturbed  surface  potential  at  the  onset 
of  saturation,  Ayf,  results  from  equating  (3.60)  and  (3.81)  with  Lqjd  =  L-L^  and 
iteratively  solving  (which  is  required  by  the  5  in  the  velocity-lateral  field  model)  the 
resulting  equation.  It  is  now  possible  to  solve  (3.84)  with  (3.85)  for  Lj  using  the 
identity  introduced  earlier.  The  accuracy  of  this  solution  is  shown  in  Figure  3.5 
which  compares  Avi/jCLqjd)  to  the  iterative  solution  with  the  discontinuity  of  the 
higher-order  derivatives  for  a  long-channel  device. 

If  the  device  is  operating  in  the  linear  region,  the  assumptions  made  for 
the  saturation  analysis  are  invalid,  but  the  resulting  equation  for  the  surface  potential 
in  the  two-dimensional  drain  region  will  be  of  the  same  form  as  (3.84): 


Vds-AVs(Lq2d) 


dy 


y  =  t,Q 


cSinh(^)+Cl2[cosh(^^]-l],  (3.86) 


where  1,,  and  C,  would  now  include  an  accounting  for  the  variation  in  inversion 
charge  density.  It  can  be  shown  that  the  two-dimensional  drain  region  is  small  in  the 


67 


1.2 
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Figure  3.5  Comparison  of  accuracy  for  the  two  solution  methods. 
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linear  region  of  operation  since   A\|;j.(Lq2J3)  =  Vpj    and   -r— 


»0   for  all 

y  =  Lq2d 


Vps  >  0.  Equation  (3.84)  is  therefore  used  for  all  bias  conditions. 

3.8  Terminal  Charge  Modeling 

The  terminal  charging/discharging  currents  result  from  the  modulation  of 
the  space  charge  as  controlled  by  the  time-dependent  terminal  biases.  These  currents 
can  be  expressed  using  the  quasi-static  approximation  as 


where  i=(G,S,D,B)  and  j=(GS,DS,BS).  Note  that  (3.87)  defines  nonreciprocal 
capacitances  and  transcapacitances  that  physically  reflect  the  charge  dynamics, 
including  short-channel  effects,  of  the  four-terminal  device. 

3.8.1  Quasi-Static  Gate  Charge 

The  gate  charge  is  defined  by  Gauss's  law  based  on  a  one-dimensional 
oxide  field,  integrated  over  the  metallurgical  channel  length: 

t-Q2D  +  t^d 

Qg  =  wc„,    J    [VGs-<i>Ms-v|/s(y)-Vp(>')l'*y>  (3.88) 

-L, 

which  is  valid  in  all  regions  of  operation.  This  integral  can  be  rewritten  as 


Qg  =  w 


CoxLlVGs-^MS-Vs(0)-v|/pWJ-Cox'      j      Av,(y)dy 


(3.89) 


where  A\(»j(y)  =  l((j(y) -^^(O)  represents  the  variation  in  the  surface  potential  with 

respect  to  the  solution  at  the  virtual  source  and  C^^'  accounts  for  the  variation  in  the 

polysilicon  depletion  potential. 

Using  the  two-dimensional  analysis  presented  in  Section  3.6,  the  integral 
can  be  evaluated  in  the  two-dimensional  source  region  resulting  in 


jAv,(y)dy  =  Cl,2[l,sinh[bj_Lj 


(3.90) 


''•'""-       J2[C„;  +  C,-Q,(0)/V,]- 

Similarly,  the  variation  in  the  potential  across  the  velocity  saturated  two- 
dimensional  drain  region  can  be  evaluated  the  using  the  two-dimensional  analysis 
presented  in  Section  3.7.  The  resulting  integral  equation  is  easily  evaluated  to  give: 


"*  L 

J      Av,(y)dy  =  Av,(Lq2b)L,  +  E^^l^^coshf^j-  l] 

+  ;i?[l,sinh(^)-L,] 


LoiD 


(3.91) 


where  L 


2(Co;  +  c,)- 
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The  variation  in  potential  in  the  quasi-two-dimensional  channel  region  is 
not  explicitly  found  during  the  calculation  of  Lq2d  and  the  channel  current.  However 
with  the  use  of  an  average  effective  mobility,  the  current  can  be  written  at  any  point 
along  the  channel  as 

'cH  =  W ^dA^«C„;  +  C,)V,  -  [Q,(0)  +  (C„;  +  C,)AvJ  }^^ ,     (3.92) 

1  -^8 T-;- 

^sat  dy 

which  can  be  integrated  from  y  =  0  to  define  the  potential  along  the  channel: 

A\j/,(y)sA-   Ia^-B-^-  (3.93) 

Qr(0)  5Ir„  ZIpHLn^n 

where  A  =  V.  -  ,--^^  -  ,,^^-^^_  and  B  =  ^  J"c  J^C.) '  '"  -'" 

to  force  the  condition  of  A\|/5(Lq2d)  =  A  -  Va--B  =  A\j(5(LQ2n) ,  B  is 
approximated  as  B  2  2AA\J/s(Lq2p)-[A\J/5(Lq2[3)]^,  which  defines  an  appropriate 
value  of  Hjff  in  (3.92).  However,  the  resulting  equation  for  the  integral  of  the 
potential  along  the  channel. 


J  Av,(y)dy  =  ALq2d  +  f  ^[(A^  -  3)3/2  _  a3]  , 


(3.94) 
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has  a  singularity  at  V^j  =  0  where  B  — >  0.  This  problem  can  be  eliminated  by  using 
a  Taylor  expansion  of  (3.94)  about  Vps  =  0  and  keeping  only  the  first  six  terms  of 
the  series: 

I 

1^02D„   .     1  ^OlD^-r,-,  .     1  ^Q 


JA..(y,dy.i::2pB.l^V.^^B3a 


l_^2Dp,4      _7_^;Q2D„5  3_^2D„6-  (3.95) 

128  a7  1536  A'  1024  aH 


3.8.2  Quasi-Static  Channel  Charge 

The  channel  charge  results  from  the  straightforward  integration  of  the 
inversion  charge  density  along  the  metallurgical  channel; 

QcH  =  W     j     Qe(y)dy.  (3.96) 

-L. 

This  integral  can  be  represented  similarly  to  the  gate  charge  with  respect  to  the 
inversion  charge  at  the  virtual  source: 


QcH  =  W 


LQ,(0)+      I      AQ,(y)dy 


(3.97) 


Consistent  with  analysis  of  the  two-dimensional  source  region,  the 
variation  in  the  inversion  charge  is  approximated  using  the  Taylor  expansion  of 
(3.75).  The  resulting  integral  is  equivalent  to  (3.90)  with  an  additional  constant  in 
front  of  the  integral: 


n 


°f  Qc(0)  °r 

jAQ,(y)dy  =  -i^  JA>)/,(y)dy.  (3,98) 


Consistent  with  the  assumption  of  constant  channel  charge  in  the  two- 
dimensional  drain  region,  the  change  in  channel  charge  in  this  region  is  simply  the 
change  in  channel  charge  at  the  boundary  multiplied  by  the  length  of  the  two- 
dimensional  region: 

J      AQ,(y)dy2(C„;  +  C,)A{fr,(LQ2D)L,.  (3.99) 

Using  the  variation  in  inversion  charge  along  the  quasi-two-dimensional 
channel  given  by  (3.59)  results  in 


Lq2D  ^32D 

J  AQ,(y)dy  =  (C„;-HC,)  j  Ax|/,(y)dy.  (3.100) 

0  0 


The  result  of  (3.95)  is  then  used  to  find  the  change  in  channel  charge  along  the 
effective  channel. 

3.8.3  Quasi-Static  Drain  Charge 

It  is  necessary  to  partition  the  channel  charge  between  the  source  and 
drain.  The  simple  scheme  presented  by  [OhSO]  has  been  chosen  to  compute  the  drain 
charge: 
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(3.101) 


-L. 


Once  again,  it  is  possible  to  represent  the  variation  in  the  drain  charge  along  the 
channel  with  respect  to  the  inversion  charge  at  the  virtual  source: 


Qd  =  w 


y  +  L3 


Qc(0)+       j      ^AQ^(y)dy 


(3.102) 


Using  the  Taylor  expansion  to  represent  the  inversion  charge  in  the  two- 
dimensional  source  region  results  in  an  expression  for  the  integral  in  (3.102)  across 
the  two-dimensional  source  region:  ;     ./  . 


y  +  L, 


J^AQ^(y)dy 


Qc(0) 


rU\    1   L? 


d.?[cosh(y-,]4 


(3.103) 


^'^"^'^=    /2[C„;.C,-Q,(0)/V,]- 

Applying  the  partitioning  scheme  to  the  two-dimensional  drain  region  results  in  an 
simple  integral  which  is  easily  evaluated  to  give 


y  +  U 


J      ^AQjy)dy  =  (C„;  +  C,)Aii>,(LQ2D)- 


2L 


(3.104) 
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Using  (3.59)  to  represent  the  variation  in  the  inversion  charge  along  the  channel,  the 
partitioning  scheme  results  in 


j  ^AQ^(y)dy  =  (C„;  +  C,) 


j  ^Av,(y)dy  +  -'  J  Av,(y)dy 


(3.105) 


The  second  integral  on  the  right-hand  side  of  (3.105)  is  evaluated  in  (3.94).  With  the 
use  of  (3.93),  the  weighted  first  integral  of  potential  can  be  approximated  as 


Y  ^Av,(y)dy  .  iAL^,„  +  ^%^[(  A^  -  B)^-]  -  ^^[(A^  -  B) 


4  Lq2D 

15  b2 


A5 


(3.106) 


As  was  the  case  in  defining  the  gate  charge,  a  Taylor  series  is  taken  about  the 
singularity  at  V^j  =  0  to  give: 


J  ^^(v)dy4^B4?^B^.l^-g5B3 

0 

5    ^2D„4         1    ^Q2Dp.5   I      21     LQ2D„f;.  (3.107) 

768  a7  256  A'  8192  A'l 


With  (3.103),  (3.104)  and  (3.107),  (3.102)  gives  Qd. 


3.8.4  Quasi-Static  Source  and  Body  Charges 

With  the  use  of  the  charge  partitioning  model,  the  source  charge  is 
expressed  as 


Lqid  +  Lj 
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Qs  =  W     J     (l  -  ^-^)Q^(y)dy 


(3.108) 


However,  this  integral  need  not  be  evaluated  since  the  source  charge  can  be 
expressed  as 

Qs  =  Qch-Qd-  (3109) 

The  body  charge  is  computed  using  the  condition  of  charge  neutrality: 

Qb  =  -(Qc  +  QcH  +  QfWL),  (3.110) 

where  Qf  is  the  fixed  charge  (areal)  density  in  the  oxide. 

3.9  Conclusions 

A  bullc  MOSFET  model  for  extremely  scaled  CMOS  which  relies  on 
structure-dependent  parameters  was  derived  to  physically  account  for  one-  and  two- 
dimensional  effects  including  charge  sharing  and  drain-induced  effects.  This  is 
uniquely  done  by  defining  a  portion  of  the  channel  where  the  two-dimensional 
effects  are  approximately  constant  and  by  using  smoothing  functions  to  define 
transitions  between  valid  physical  approximations.  The  quasi-two-dimensional 
channel  length,  Lq2d,  is  bias-dependent  allowing  reliable  prediction  of  current  and 
other  device  characteristics.  UFET  is  C^ -continuous  and  process-based,  which 

makes  it  an  ideal  candidate  for  process  integration,  statistical/sensitivity  studies  and 
circuit  simulation  in  extremely  scaled  CMOS  technologies. 


CHAPTER  4 

UFET  IMPLEMENTATION,  PARAMETER  EVALUATION  AND 

VERIFICATION 


4.1   Introduction 

The  new  UFET  model  formalism  presented  in  Chapter  3  along  with  the 
LDD  analysis  and  local-field  impact-ionization  model  of  Chapter  2  has  been 
implemented  in  both  MMSPICE-4.0  and  a  simulator-independent  API  (Application 
Program  Interface).  After  the  model  implementation  is  discussed,  the  intrinsic  model 
parameters  introduced  in  the  model  formalism  and  the  additional  models  for  the 
extrinsic  circuit  elements  are  presented.  Proper  evaluation  of  these  model  parameters 
is  necessary  in  order  to  reliably  and  efficiently  use  UFET.  A  method  of  obtaining  the 
process-based  model  parameters,  which  relies  on  knowledge  of  the  device  physics  in 
all  regions  of  operation  and  of  device  structure,  is  presented  and  highlights  the 
usefulness  of  process-based  modeling  where  the  model  parameters  correlate  with  the 
actual  process.  This  methodology  is  demonstrated  with  a  typical  retrograde  device 
simulated  in  MEDICI  [Tec92]  across  channel  length  and  bias.  After  a  comparison  to 
the  simulated  DC  and  AC  data,  UFET  is  applied  to  devices  in  two  deep-submicron 
technologies  to  assess  the  model  against  measured  data. 

4.2  UFET  Implementation 

The  network  representation  of  UFET  is  shown  in  Figure  4.1.  The 
implementation  of  the  model  involved  developing  a  subroutine  (SSRMOD)  to  effect 
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the  model  formalism  as  described  in  Chapter  3.  The  structure  of  the  intrinsic  model 
routine  is  shown  in  Figure  4.2.  Since  UFET  does  not  use  a  regional  approach,  as  in 
[Gre93a],  the  boundaries  for  moderate  inversion  are  not  needed,  and  the  channel 
current  and  terminal  charges  are  calculated  using  single  functions  which  improves 
the  continuity  and  convergence  properties. 

4.2.1  MMSPICE 

The  new  model  formalism  presented  in  Chapter  3  (UFET)  has  been 
implemented  in  MMSPICE  [Jeo90,  Gre93a]  version  4.0,  an  enhanced  version  of 
SPICE2  [Nag75]  which  contains  a  physical  charge-based  BJT/HBT  model  [Hon94] 
in  addition  to  the  standard  device  models.  The  model  is  implemented  as  a  series  of 
subroutines  (see  Figure  4.3)  which  calculate  the  equivalent  currents  and 
conductances  of  a  MOSFET.  QBFET,  which  is  similar  to  the  subroutine  MOSFET  in 
the  original  SPICE2  code,  calls  routines  that  analyze  the  intrinsic  (QBID)  and 
extrinsic  (QBEC)  regions  of  a  device.  QBID  in  turn  calls  SSRMOD  repeatedly  to  find 
the  currents  (I,;^  and  I^j),  the  conductances  and  transconductances  of  the  currents 

and  the  capacitances  and  transcapacitances  associated  with  the  intrinsic-region 
charges  (QG,  QB,  QS  and  QD),  using  a  finite-difference  method.  The  extrinsic 
charge  storage  is  then  accounted  for  in  QBEC. 

4.2.2  Application  Program  Interface 

In  the  course  of  evaluating  the  UFET  model,  several  deficiencies  (e.g., 
lack  of  AC  modeling  and  transient  simulation  problems)  of  the  underlying 
implementation  in  MMSPICE  became  apparent.  It  was  determined  that  the  best  of 
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Figure  4.2  Flowchart  of  model  algorithm  (SSRMOD). 
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course  of  action  was  to  implement  UFET  in  a  standard  (simulator  independent)  API 
which  allows  different  simulator  platforms  to  use  identical  model  source  code.  Using 
this  Type-I  interface  [New95],  UFET  was  glued  to  SPICE3e2  [Qua92]. 

The  Type-I  interface  is  designed  to  capture  the  critical  elements  of  a 
MOSFET  model  in  a  single  copy  of  the  source  code  which  can  be  easily  distributed 
and  installed  on  any  simulator  platform  as  shown  in  Figure  4.4.  This  is  accomplished 
(see  Figure  4.5)  by  defining  a  set  of  data  structures  to  hold  model  and  instance  data 
and  using  a  predefined  subroutine  name  for  the  various  calculations  necessary  in  a 
MOSFET  model. 

Since  the  interface  is  restricted  to  a  four-terminal  MOSFET  with  a  pre- 
defined matrix  stamp,  several  effects  can  not  be  included.  For  instance,  the  parasitic 
source/drain  resistances  must  be  bias-independent,  and  there  is  no  provision  made  for 
modeling  self-heating  or  non-quasi  static  effects.  Within  these  limitations,  the 
interface  still  provides  advantages  over  the  standard  method  of  hard-wiring  a  model 
into  a  simulator.  The  interface  allows  the  distribution  of  simulator-independent 
model  code  as  either  source  or  object  code  which  guarantees  that  the  same  model  is 
being  used  across  platforms.  This  also  simplifies  the  work  needed  to  upgrade  the 
models  as  research  continues. 

4.3  Instance  and  Model  Parameters 

The  model  is  based  on  eight  instance  parameters  (see  Table  4.1)  and  31 
model  parameters  shown  in  Table  4.2.  The  instance  parameters  vary  among 
transistors  in  a  technology  or  circuit  and  are  the  gate  length  (L)  and  width  (W),  the 
areas/perimeters  of  the  source  (AS,  PS)  and  drain  (AD,  PD),  the  device  temperature 
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Figure  4.4  Diagram  of  the  Type-I  standard  interface  to  MOSFET  models. 
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Figure  4.5  Design  of  the  Type-I  standard  interface  to  MOSFET  models. 
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Table  4.1   MOSFET  Instance  Parameters 


Name 

Description 

Units 

Default 

Typical 
Value 

ND 

Drain  Node 

- 

- 

- 

NG 

Gate  Node 

- 

- 

NS 

Source  Node 

- 

- 

- 

NB 

Body  Node 

- 

- 

- 

MNAME 

Model  Name 

- 

- 

- 

L 

Gate  length 

m 

l.Oe-6 

0.1 e-6 

W 

Gate  width 

m 

l.Oe-5 

I.Oe-5 

AD 

Drain  area 

m2 

0 

I.Oe-U 

PD 

Drain  perimeter 

m 

0 

AS 

Source  area 

m^ 

0 

l.Oe-lI 

PS 

Source  perimeter 

m 

0 

M 

Multiple  device  option 

- 

1 

1 

TEMP 

Device 
temperature 

"C 

Circuit 
Temperature 

Off 

DC  analysis  with  device  off  flag 

- 

- 

- 

IC 

Device  initial  conditions 

V 

- 

- 

Note:  The  parameters  PD,  PS,  M,  and  TEMP  are  not  implemented  in  MMSPICE, 
but  have  been  implemented  in  the  API. 
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Table  4.2  MOSFET  Model  Parameters 


Name 

Description 

Units 

Default 

Typical 
Values 

DL 

Channel-length  reduction 

m 

0.0 

0.1 e-6 

DW 

Channel-width  reduction 

0.0 

0.1 e-6 

NB 

Surface  channel  doping  density 

cm 

1.0el7 

1.0el7 

NBH 

High  side  doping  density 

cm-3 

l.OelS 

1.0el8 

XLH 

Position  of  high/low  junction  in 
retrograde  device 

m 

1.0 

0.05e-6 

XJ 

Source/Drain  junction  depth 

m 

0.1 e-6 

0.1 e-6 

TOX 

Gate  oxide  thickness 

m 

5e-9 

5e-9 

TPG 

Type  of  gate  material: 

-Hi)  opposite  to  body 

-1)  same  as  body 

0)  Aluminum 

- 

+  \ 

+  1 

NPOLY 

Polysilicon-gate  doping  density 

cm-3 

le25 

5el9 

NQF 

Oxide  fixed  charge 

cm"^ 

0.0 

l.Oe-10 

WKF 

Gate  work  function  difference 

V 

Calculated 

- 

LLDD 

LDD  region  length  (0  for  no  LDD) 

m 

0.0 

0.15e-6 

NLDD 

LDD  doping  density 

cm-3 

0.0 

lel8 

RLDD 

Specific  LDD  parasitic  resistance 

n-m 

0.0 

500e-6 

UO 

2^ro-field  mobility 

cm^/(Vs) 

Calculated 

700  (NMOS) 
400  (PMOS) 

THETA 

Linear  mobility  degradation 
coefficient 

cntTV 

l.Oe-6 

I.Oe-6 

THETA2 

Quadratic  mobility  degradation 
coefficient 

cm^A'^ 

I.Oe-12 

l.Oe- 12 

FB 

Vps-averaging  factor  for  mobility 
degradation 

- 

1.0 

1.0 

VSAT 

Carrier  saturated  drift  velocity 

cm/s 

l.Oe? 

l.Oe? 

Table  4.2  —  continued 


Name 

Description 

Units 

Default 

Typical 
Values 

Al 

Impact-ionization  parameter  Uq 

cm-' 

0.0 

2.45e6  (NMOS) 
3.7e6  (PMOS) 

Bl 

Impact-ionization  parameter  Pq 

V/cm 

0.0 

l.92e6{NMOS) 
8.0e6  (PMOS) 

QM 

Quantum  mechanical  effects 
modulation  factor 

- 

1.0 

1.0 

FCS 

Short-channel 
charge-sharing  modulation  factor 

- 

0.0 

0.1 

NO 

Peak  doping  in  RSCE/halo 

cm-3 

0.0 

lel8 

YO 

Characteristic  length  of 
RSCE/halo 

m 

0.0 

0.05e-6 

CGDSO 

Gate-(source  -H  drain)  overlap 
capacitance 

F/m 

Calculated 

l.Oe-10 

CGBO 

Gate-body  overlap  capacitance 

F 

0.0 

l.Oe-0 

RDS 

Total  specific  extrinsic  resistance 

n-m 

0.0 

600e-6 

JS 

Parasitic  diode  saturation  current 
density 

A/rn^ 

0.0 

I.Oe-9 

N 

Parasitic  diode  ideality  factor 

- 

1.0 

1.0 

CJO 

Zero-bias  parasitic  diode  junction 
capacitance 

F/m^ 

0.0 

l.Oe-3 

MJ 

Parasitic  diode  junction  grading 
factor 

- 

0.5 

0.5 

PJ 

Parasitic  diode  junction  potential 

V 

0.8 

0.8 

Note  on  differences  between  MMSPICE  and  API  implementations:  The 
parameters  CGDSO  and  RDS  have  replaced  CGDO,  CGSO,  RD  and  RS 

respectively.  In  MMSPICE,  the  extrinsic  overlap  capacitance  and  resistance  should 
be  divided  between  the  respective  source  and  drain  components.  The  parameters 
LLDD  and  NLDD  have  replaced  LLD  and  NLD  as  used  in  MMSPICE.  Also,  the 
parameter  ULD  has  been  replaced  by  RLDD  which  is  the  quantity  actually  measured/ 
evaluated. 
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(TEMP)  and  the  multiple  finger  option  (M).  The  31  model  parameters  are  assumed 
to  be  constant  for  devices  in  a  given  technology  and  describe  the  intrinsic  and 
extrinsic  portions  of  the  device.  Of  the  31  model  parameters,  20  are  physically 
significant  and  can  be  extracted  from  the  device  cross-section  or  doping  profiles  of 
the  device. 

4.3.1  Intrinsic  Model  Parameters 

The  intrinsic  portion  of  the  device  is  described  by  the  23  parameters  which 
define  the  structure,  (approximate)  doping  profiles  and  doping  dependent 
parameters.  Of  these  23  parameters,  only  seven  can  not  be  estimated  from 
information  about  the  process  but  must  be  input  by  the  user.  The  carrier  saturation 

velocity  (VSAT)  is  a  material  constant  with  a  typical  value  of  -  8x10  cm/s .  The 
impact-ionization  model  relies  on  two  empirical  constants  (AI  and  BI)  in  order  to 
calculate  the  impact-ionization  rate  in  the  channel  and  LDD.  The  transverse-field 
dependence  of  the  carrier  mobility  in  the  channel  is  modeled  by  a  linear  (THETA) 
and  quadratic  (THETA2)  dependence  on  the  vertical  field,  although  a  recent  paper 
presented  a  derivation  of  these  two  parameters  from  basic  scattering  processes 
[Rei96].  The  carrier  mobility  is  also  a  function  of  the  potential  along  the  channel,  and 
an  approximate  integral  is  used  with  a  fitting  parameter  (FB)  to  determine  the 
channel  current. 

4.3.2  Extrinsic  Model  Parameters 

The  remaining  8  parameters  describe  the  extrinsic  portion  of  the  device. 
In   scaled  devices,   it  is  a  necessity  to  accurately   model   the  device  extrinsics 
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(including  interconnects  which  are  not  part  of  UFET)  in  order  to  accurately  model 
transient  and  AC  device  behavior. 

The  parasitic  junction  diodes  each  have  five  model  parameters  with  the 
junctions  assumed  to  be  symmetric.  The  leakage  current  through  the  parasitic  source 
diode  is  [Ant88]: 


(4.1) 


where  JS  is  the  junction  leakage  current  and  MJ  is  the  junction  grading  coefficient. 
The  parasitic  junction  capacitances  are  assumed  to  be  predominantly  due  to  the  area 
component.  The  depletion  charge  associated  with  the  source  junction  is  [McA93] 


BS 

=  AS  JS 

1 

^ 

/ 

Qbs  =  AS^,(.-^)  ,  ,  (4.2) 


with  CJO  being  the  zero-bias  depletion  capacitance  associated  with  the  junction,  PJ 

Inl  I  j-pA(R*PJ-Vbs)i 

IS     the     built-in    junction    potential     and     Vbs  =  R*PJ-     '        '- 

A 

continuously  models  the  voltage  dependence.  The  constants  R  and  A  are  generally 

extracted  under  forward-bias  conditions  but  have  been  hard-wired  to  reasonable 

values. 

The  overlap  capacitors  are  modeled  as  simple  parallel  plate  capacitors 

(CGDSO  and  CGBO)  and  contribute  to  the  charging  currents  shown  in  Figure  4.1. 

In  addition  to  the  voltage  drops  across  the  LDS/LDD  regions  of  the  device,  the 

additional  extrinsic  resistances,  e.g.,  contact  resistance,  are  accounted  for  with  the 


WW 
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model  parameter  RDS  (=  Rj  +  Rq).  In  order  to  improve  the  computation  efficiency 
of  tlie  implementation,  the  ohmic  portions  of  the  LDS/LDD  drops  (Rlds  and  Rldd) 
are  added  to  Rj  and  R^  respectively  as  shown  in  Figure  4. 1 . 

4.3.3  Model  Parameter  Temperature  Dependencies 

Since  the  model  parameters  are  physical,  the  temperature  dependencies  of 
the  parameters  are  implemented  using  published  models  [Wor98]  in  a 
straightforward  manner.  An  example  of  this  is  the  model  for  the  intrinsic  carrier 
density 

n^  =  l.45xl0"*(3^]''e<'^E,  ^4  3^ 


-4      T^ 
1.16-7.02x10 


u           AC            1115                                  T+  1108       ,,  . 
where     AEq  =  ^^^ 2kT '        ^'"^  temperature- 

dependent  value  of  Uj  leads  to  an  implicit  temperature-dependent  threshold  voltage 

through  the  definition  of  (|),^  =  V.ln^ 


"■      ■    'n. 


4.4  General  Approach  to  Parameter  Evaluation 

Proper  evaluation  of  the  model  parameters  is  necessary  in  order  to  reliably 
and  efficiently  use  UFET.  Since  most  of  the  parameters  are  based  either  on  the  device 
structure  or  the  underlying  physics,  the  resulting  parameter  evaluation  methodology 
is  relatively  straightforward  and  does  not  rely  on  a  standard  optimizer.  The  key  to 
this  approach  is  to  recognize  that  a  limited  subset  of  the  model  parameters  controls 
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the  DC  or  AC  terminal  characteristics  in  a  given  region  of  operation.  The  majority  of 
the  model  parameters  can  be  evaluated  from  the  DC  measured  characteristics. 

Table  4.3  Technological  MOSFET  Parameters  Used  for  Verification 


Parameter 

Value 

TOX 

5.5  nm 

DL 

0.05  nm 

XJ 

0.2  nm 

L 

0.35  \lm 

W 

20|im 

To  start  with,  an  initial  set  of  parameters  is  assumed,  from  processing 
(technology)  or  2-D  process  simulation,  e.g.,  with  SUPREM  [Tec89].  Since  UFET  is 
physical  with  several  structural  parameters,  this  initial  evaluation  is  straightforward 
and  indeed  can  give  very  good  estimates  for  the  parameter  values.  The  (assumed) 
initial  values  for  the  parameters  however  need  not  be  known  with  a  great  degree  of 
certainty.  Once  the  initial  values  are  obtained,  the  calibration  methodology  can  be 
applied  in  a  systematic  manner  using  a  short-  and  a  long-channel  device.  This 
methodology  is  demonstrated  with  a  typical  retrograde  device  simulated  in  MEDICI 
[Tec92]  across  channel  length  and  bias.  The  retrograde  MOSFET  used  for  the 
comparison  is  based  on  a  device  presented  at  a  recent  lEDM  [Ven95]  with  the 
structural  and  doping  parameters  shown  in  Table  4.3  and  the  channel  profile  shown 
in  Figure  4.6. 
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Figure  4.6  Retrograde  channel  doping  profile  used  in  MEDICI. 
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4.5  Parameter  Evaluation  for  Retrograde  Bulk  MOSFETs 

The  model  parameters,  along  with  their  description  are  listed  in  Table  4.2. 
Also  listed  are  their  typical  values  for  current  state-of-the-art  bulk  technologies. 
Now,  given  the  technology,  TOX  is  usually  well  known,  and  hence  it  can  be  directly 
specified  in  the  input  model  card.  The  critical  parameters  in  the  table  that  usually 
have  to  be  evaluated  are  XLH,  NB,  NBH,  NPOLY,  DL,  RDS,  UO,  THETA, 
THETA2  and  VSAT.  If  an  LDS/LDD  is  present  in  the  structure,  LLDD  and  NLDD 
must  also  be  estimated.  Some  of  the  other  parameters  (e.g.,  FB)  have  only  a 
secondary  effect  on  the  device  characteristics.  Hence,  parameter  evaluation  is 
restricted  to  only  a  few  structural  and  physical  parameters.  The  evaluation  process 
will  be  systematically  detailed  by  considering  different  regions  of  operation  where 
these  key  parameters  have  characteristic  signatures. 

4.5.1   Results  from  Assumed  Initial  Parameter  Values 

Based  on  the  channel  profile,  the  initial  parameter  values  can  be  used  to 
start  the  parameter  evaluation  process.  The  value  of  NBH  is  assumed  to  be  equal  to 
the  peak  doping  in  the  channel,  and  NB  is  set  equal  to  a  value  between  the  surface 
doping  concentration  and  NBH.  The  position  of  the  high-low  junction  is  then 
estimated  to  be  one-half  of  the  position  of  the  peak  doping  in  the  channel.  Using  the 
standard  mobility  curves  then  provides  an  estimate  of  UO.  The  other  parameter 
values  are  then  set  to  the  UFET  default  values.  The  parameter  values  used  for  this 
first  evaluation  are  shown  in  Table  4.4  in  the  column  labeled  Level  =  0.  Predictions 
of  current-voltage  characteristics  resulting  from  the  use  of  this  initial  parameter  set 
are  shown  in  Figures  4.7-4.10.  It  should  be  noted  that  the  large  error  for  the  short- 
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Figure  4.7  Subthreshold  plot  with  Vgs  =  OV  using  Level  =  0  model  parameters. 
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Figure  4.8  Subthreshold  plot  with  Vgj  ^  OV  using  Level  =  0  model  parameters. 
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Figure  4.9  Threshold  plot  with  V^s  =  OV  using  Level  =  0  model  parameters. 
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Figure  4. 1 0  Linear-Saturation  plot  with  Vbs  =  OV  using  Level  =  0  model  parameters. 


channel  device  in  the  threshold  plot  of  Figure  4.9  is  due  to  the  default  ideal  contact 
in  MEDICI.  Although  a  better  initial  estimate  of  the  parasitic  resistances  would  have 
been  zero  in  this  case,  the  first  stages  of  the  parameter  evaluation  are  done  with  the 
UFET  default  values  in  order  to  demonstrate  the  initial  model  parameter  values  need 
not  be  very  accurate  to  calibrate  the  model. 

4.5.2  Calibration  Based  on  Subthreshold  Characteristics 

In  the  subthreshold  region,  the  current  conduction  mechanism  in  a 
MOSFET  is  predominantly  by  carrier  diffusion  [SzeSl],  and  can  hence  be  modeled 
in  a  rather  straightforward  manner.  The  current  in  this  region  can  then  be  used  to 
obtain  information  on  the  position  of  the  high/low  junction  XLH,  the  (uniform)  body 
doping  NB  and  NBH  (for  Vgs  *  0) ,  and  the  lateral  source/drain  diffusion  length  DL. 

Now,  XLH,  DL,  NB  and  NBH  are  calibrated  for  a  given  device  at  low 
drain  bias.  The  subthreshold  slope  is  given  approximately  as 

Ss6ofl+-il  (4.4) 


where    Cj  s     J^       is  the   depletion   capacitance   and   C^,^  =       ""      is   the   oxide 
XLH  1  OX 

capacitance.   With   the   initial   estimates   of  the   doping  concentrations,   DL   and 

knowledge  of  TOX  from  the  technology,  XLH  can  then  be  tuned  to  obtain  the 

subthreshold  slope.  NB  and  NBH  usually  have  only  a  minor  effect  on  subthreshold 

slope.  Using  this  value  of  XLH  gives  in  the  results  of  Figure  4.11  which  use  the 

parameters  shown  as  Level  =  1  in  Table  4.4. 


W  =  20|^m 


Figure  4. 1 1   Subthreshold  plot  with  Vgs  =  OV  using  Level  =  1  model  parameters. 


100 

The  current  in  this  region,  which  depends  on  the  gate-controlled  body 
charge  (Qg  <=<:  NB*XLH) ,  can  now  be  used  to  calibrate  NB.  If  necessary,  XLH  may 

also  be  fine-tuned  to  fit  both  the  current  and  the  slope.  Note  that  usually  the  lateral 
source/drain  diffusion  length  DL  is  a  constant  for  all  channel  lengths  in  a  given  wafer 
even  when  RSCE  is  significant.  To  ensure  consistency,  the  calibration  procedure 
may  be  applied  to  a  range  of  channel  lengths  to  define  an  "optimized"  DL  for  the 
given  technology  based  on  the  current  as  well  as  the  slope  at  low  and  high  Vqj.  The 
parameters  as  this  stage  of  the  evaluation  process  are  shown  in  Table  4.4  as  Level  = 
2,  and  the  corresponding  plot  is  shown  in  Figure  4.12. 

If  data  is  available  for  Vgj  ?^0,  the  doping  concentration  on  the  high-side 

of  the  junction  can  also  be  determined  from  the  subthreshold  data.  The  potential  drop 
across  the  high  side  of  the  junction  results  in  a  reduction  of  the  body  bias  effect 
controlled  by  the  doping  concentration.  By  tuning  the  doping  concentration  in  this 
region  of  the  device  across  several  Vgs's,  an  "optimized"  value  of  NBH  can  be 
found.  If  the  optimized  value  of  NBH  is  significantly  different  from  the  initial 
estimate,  the  value  of  NB  may  also  need  to  be  slightly  tuned  in  order  to  improve  the 
fit  at  low  body  bias.  This  was  true  in  this  case  due  to  the  large  value  of  NBH  initially 
assumed,  and  this  is  reflected  in  the  Level  =  3  parameters  of  Table  4.4.  The  results 
of  the  body  bias  simulations  are  shown  in  Figure  4.13 

At  this  stage  of  the  tuning  process,  XLH,  DL,  NB  and  NBH  are  known. 
The  remaining  parameters  are  now  extracted  in  strong  inversion. 
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Figure  4.12  Subthreshold  plot  with  Vgs  =  OV  using  Level  =  2  model  parameters. 
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Figure  4.13  Subthreshold  plot  with  Vgs  -^  OV  using  Level  =  3  model  parameters. 
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4.5.3  Calibration  Based  on  Strong-Inversion  Characteristics 

In  strong  inversion,  the  model  parameters  which  can  be  evaluated  are  UO, 
THETA,  THETA2,  RDS  and  VSAT.  If  an  LDS/LDD  is  present,  LLDD  and  NLDD 

can  also  be  evaluated  here.  For  all  structures,  it  is  assumed  that  the  parasitic  source 
and  drain  resistances  are  equal  (and  so  are  Lmj  and  L^d^)  due  to  symmetry. 

At  low  Vqs,  the  linear  region  characteristics  of  the  MOSFET  are  similar 
to  that  of  a  resistor  with  its  total  ON  resistance  (Rqn)  g'ven  approximately  by 


^  ^^  ^  RDS     ^;^ff , .  ,, 

ON-ij^g-    w       WC„,(VQg-V.J,)^(UO,  THETA,  THETA2)'        '' 


As  can  be  seen  from  (4.3),  Rqn  's  dependent  on  UO,  THETA,  THETA2  and  RDS. 
To  isolate  the  effects  of  RDS  from  those  of  UO/THETA/THETA2,  it  can  be  shown 

2 

from  (4.5)  that  g^/Ios  '^  nearly  independent  of  RDS.  This  allows  us  to  first  evaluate 

UO,  THETA  and  THETA2  to  fit  the  g^/Ios  characteristics,  independent  of  RDS. 
Once  these  three  parameters  are  tuned,  RDS  can  be  explicitly  calibrated  from  (4.5). 
If  an  LDS/LDD  is  present  in  the  structure,  LLDD  is  set  equal  to  the  spacer  width  as 
defined  by  technology,  and  NLDD  is  set  to  roughly  correspond  to  the  LDD  dose. 
Typically  the  LDD  doping  is  in  the  mid-  10'^  cm"'  range.  When  an  LDS/LDD  is 
present  the  resistance  RDS  in  (4.5)  would  include  implicitly  the  ohmic  component  of 
the  series  resistance  (2*RLDD).  However,  the  RDS  in  the  model  card  is  the  source/ 
drain  parasitic  resistance  without  RLDD.  Hence,  2*RLDD  must  be  subtracted  off 
from  the  evaluated  (total)  RDS.  Note  that  if,  as  in  some  technologies,  the  LDS/LDD 


104 

is  more  of  a  source/drain  extension,  then  the  effect  of  the  LDS/LDD  is  to  simply 
increase  the  total  resistance  RDS  and  the  non-ohmic  drops  in  the  LDD  will  not  be 
important.  For  such  cases,  it  is  recommended  that  the  LDS/LDD  feature  be  turned  off 
in  the  model  to  speed  up  run  time.  The  extraction  technique  for  RDS  described  above 
should  then  be  carried  out  as  though  there  were  no  LDS/LDD.  After  evaluating  these 
parameters  (Level  =  4  in  Table  4.4),  the  resulting  threshold  plot  is  Figure  4,14. 

At  high  Vqs,  and  relatively  large  Vqj.  the  current  saturates  in  short- 
channel  MOSFETs  due  to  carrier  velocity  saturation.  From  the  model  formalism,  the 
current  for  a  short-channel  device  in  strong  inversion  at  high  Vpj  is 
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From  (4.6)  and  the  output  Ids'^ds  characteristics,  VSAT  can  be  evaluated  at  V^s 
close  to  VQS(sai)  the  saturation  drain  voltage  at  moderately  high  Wqs-  (^^e  initial 

guess  for  VSAT  can  be  set  to  6x10  cm/s  for  NMOS  although  1x10  cm/s  was  used 
in  order  to  be  consistent  with  MEDICL)  The  final  parameter  set  (Level  =  5  in 
Table  4.4)  results  in  the  linear-saturation  characteristics  of  Figure  4.15. 

At  high  Vps,  the  current  saturates  in  long-channel  MOSFETs  due  to 
carrier  pinch-off  at  the  drain  end.  From  the  model  formalism,  the  current  for  a  long- 
channel  device  in  strong  inversion  at  high  Vps  is 
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Figure  4.14  Threshold  plot  with  Vgs  =  OV  using  Level  =  4  model  parameters. 
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Figure  4. 1 5  Linear-saturation  plot  with  Vbs  =  OV  using  Level  =  5  model  parameters. 
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From  (4.7)  and  the  long-channel  output  Ids'^DS  characteristics,  the  values  of  XLH 
and  NBH  which  both  affect  Cj  can  be  further  tuned  from  the  values  obtained  in  the 

subthreshold  region  of  operation. 

■'   '[  '■''  ' 

4.5.4  Short-Channel  Effects  Parameter 

In  scaled  devices,  the  short-channel  effects  significantly  affect  the  device 
output  characteristics.  A  single  parameter  is  available  to  tune  the  average  short- 
channel  effects  (FCS).  FCS  controls  the  Vqs  =  0  short-channel  effects  and  should  be 
tuned  in  the  low-drain  bias  subthreshold  region  for  a  minimum  channel  length 
device.  This  parameter  accounts  for  the  non-parabolic  potential  which  was  assumed 
in  deriving  an  expression  for  the  two-dimensional  effects. 

4.5.5  Reverse  Short-Channel  Effect  Parameters 

In  scaled  devices,  it  is  also  necessary  to  account  for  nonuniform  lateral 
doping  profiles  due  to  hales  and  also  the  Reverse  Short-Channel  Effect  (RSCE).  This 
is  accounted  for  by  using  a  lateral  doping  transform  which  assumes  quasi-Gaussian 
profiles  (with  peak  doping  of  Nj,  and  a  characteristic  length  y,,)  near  the  source/drain 
junctions  modify  the  surface  doping  concentration  as  a  function  of  the  channel  length 
[Bru95].  The  two  parameters  associated  with  this  effect  should  be  tuned  to  the 
threshold  voltages  of  two  short-channel  devices  showing  RSCE. 
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4.6  Application  to  Simulated  0.35^m  Technology 

4.6.1  0.35nm  Technology 

With  the  UFET  MMSPICE/API  implementations,  verification  is  provided 
by  comparing  subthreshold,  threshold  and  linear-saturation  characteristics  over 
different  channel  lengths  and  substrate  biases  with  numerical  simulations  using 
MEDICI.  For  the  verification,  the  parameters  evaluated  in  Section  4.5  are  used,  and 
it  should  be  stressed  that  no  attempt  was  made  to  find  the  optimal  values  for  the 
parameters.  It  is  also  important  to  note  that  the  same  parameters  are  used  for  all 
simulation  results  shown  in  this  section  without  "tweaking"  to  give  the  best  results 
for  each  individual  figure. 

4.6.2  DC  Verification 

The  low- Vqs  behavior  is  exhibited  in  the  threshold  plot  of  Figure  4. 1 6  and 
the  subthreshold  plot  shown  in  Figure  4.17  with  no  substrate  bias.  With  parameters 
tuned  from  the  0.3nm  and  10|im  devices  in  Figure  4.5,  the  model  accurately  predicts 
the  channel  length  and  Vqs  dependence  of  the  channel  current  for  channel  lengths 
ranging  from  long  to  shorter  than  the  target  length  (0.3nm)  for  which  the  technology 
is  designed.  The  Vgs  dependence  of  the  model  is  shown  in  Figures  4. 1 8  and  4. 1 9  for 
the  nominal  channel  length  of  the  technology.  The  channel  current  is  accurately 
modeled  in  all  regions  of  operation  by  accounting  for  the  increased  doping  in  the  p* 
region  along  the  channel,  and  UFET  also  shows  a  continuous  transconductance  result 
in  Figure  4.20.  The  discrepancies  in  the  V3S  plots  can  be  attributed  to  the  different 
mobility  models  used  in  MEDICI  and  UFET,  and  the  transconductance  errors  are 
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Figure  4.16  Threshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
for  Vgs  =  OV  showing  channel  length  dependence. 
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Figure  4.17  Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
for  Vgs  =  OV  showing  channel  length  dependence. 
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Figure  4.18  Threshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
for  Lgff  =  0.3nm  showing  ¥33  dependence. 
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Figure  4.19  Subthreshold  of  UFET-  and  MEDlCI-simulated  data  plot  for  L^ff  = 
0.3|xm  showing  Vgs  dependence. 


113 


L.„  =  0.3[im      W  =  20|im 


2.0 


1.5 


0.5 


• 

i   *  \ 

r  *  • 

'  '  '  1  1  ' 

■ 

•Vbs=   ov 

•>* 

■Vbs  =  -1V 

■rK-* 

♦V„,  =  -2V 

/   /^  * 

. 

Vos  =  0.1V 

1          ^^^ 

• 

1          v^ 

■ 

\8\ 

:i 

' 

- 

^ 

V. 

■ 

\»\. 

♦ 

\^v    ; 

\^bn^ 

<^^ 

11 

{»»^»--^ 

,     ,     ,     1     ,     ,     ,     , 

0.0  0.5 


1.0 
V 


1.5  2.0  2.5 


GS 


(V) 


Figure  4.20  Transconductance  plots  of  UFET-(lines)  and  MEDICI-(points)simulated 
data  for  L^ff  =  0.3|iin  showing  Vbs  dependence. 
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what  result  from  the  parameter  evaluation.  In  order  to  better  model  the 
transconductance,  the  mobility  parameters  (UO,  THETAl  and  THETA2)  should  be 
evaluated  from  a  transconductance  plot.  The  Vq^  dependence  of  the  model  is  shown 
in  Figure  4.21.  The  channel  current  is  well  modeled  for  both  the  long-  and  short- 
channel  cases  for  low  and  high  drain  biases.  The  long-channel  device  shows  no  drain- 
induced  effects  as  expected,  and  the  short-channel  device  shows  significant  effects 
from  DIBL  in  the  subthreshold  region  of  operation.  The  error  in  the  subthreshold 
slope  for  the  short-channel  device  at  high  drain  bias  can  be  attributed  to  deviations 
from  the  assumed  parabolic  potential  used  to  define  the  two-dimensional  effects. 

The  long-channel  linear-saturation  characteristics  shown  in  Figure  4.22 
clearly  shows  the  pinch-off  behavior  of  the  channel  current  of  the  device  for  all  gate 
biases.  The  model  gives  good  fits  to  the  MEDICI  data  since  the  non-fundamental 
parameters  such  as  the  saturated  drift  velocity  and  the  mobility  dependence  on  lateral 
field  only  play  a  minor  part  in  determining  the  output  characteristics  of  this  long- 
channel  device.  The  significant  error  the  long-channel  device  is  due  to  the 
assumption  of  a  constant  depletion  capacitance  along  the  channel  to  account  for  the 
depletion  charge  variation.  In  a  long-channel  device  controlled  by  pinch-off,  this 
assumption  tends  to  under-predict  the  saturation  current,  and  the  magnitude  of  the 
error  is  dependent  on  well  the  actual  doping  profile  tends  to  keep  the  depletion 
capacitance  constant  along  the  channel.  For  the  nominal  channel  length  device  shown 
in  Figure  4.23,  the  effect  of  channel  length  modulation  due  to  the  applied  drain  bias 


is  clearly  evidenced  by  the  increase  in  the  drain  conductance 
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Figure  4.21  Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
for  Vbs=OV  showing  DIBL. 
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Figure  4.22  Linear-saturation  plots  of  UFET-{Iines)  and  MEDICI-(points)simulated 
data  for  L^ff  =  9.95|xm  with  Vgs  =  OV. 
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Figure  4.23  Linear-saturation  plots  of  UFET-(lines)  and  MEDICI-(points)simulated 
data  for  L^ff  =  0.3nm  with  V35  =  OV. 
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compared  to  the  long-channel  case.  There  is  a  good  fit  for  the  channel  current  with 
some  error  in  the  drain  conductance,  for  which  there  is  no  parameter  in  the  model  as 
shown  in  Figure  4.24.  The  likely  cause  for  the  error  in  the  drain  conductance  is  the 
DL  assumed  for  the  devices  in  this  technology.  Since  the  technology  was  simulated 
with  MEDICI,  the  DL  was  assumed  to  represent  the  difference  between  the  physical 
gate  length  and  the  metallurgical  channel  length,  but  a  recent  paper  [Niu97]  indicates 
that  the  effective  (i.e.,  gate-controlled)  channel  length  is  larger  than  the  metallurgical 
channel  length  suggesting  that  DL  should  be  smaller  than  the  value  used.  This  will 
decrease  the  drain  conductance  while  necessitating  an  increase  in  VSAT  which  tends 
to  be  different  from  the  MEDICI  value  due  to  the  differing  mobility  models  used. 
This  dependence  of  drain  conductance  on  DL  implies  the  possibility  that  DL  should 
be  tuned  to  best  fit  the  drain  conductance  if  this  does  not  compromise  the  model 
accuracy  in  other  regions  of  operation. 

4.6.3  Charge  Model  Verification 

The  DC  portion  of  the  model  has  been  validated  against  MEDICI 
simulations,  but  the  new  charge  modeling  must  now  be  verified.  The  charge  model 
is  qualitatively  verified  by  looking  at  the  intrinsic  capacitances  and  trans- 
capacitances  in  both  long  and  short-channel  devices.  These  intrinsic  capacitances  are 
defined  by  the  following  relationships. 
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Figure  4.24  Drain  conductance  plots  of  UFET-(lines)  and  MEDICI- 
(points)simulated  data  for  Lj.ff  =  0.3nm. 
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The  UFET  intrinsic  gate  capacitances  for  a  long-channel  device  are  shown 
in  Figure  4.25.  Notice  the  smooth  transitions  between  the  different  regions  of 
operation,  i.e.,  accumulation/depletion,  subthreshold/threshold  and  linear/saturation, 
that  result  from  the  C_  continuity  of  the  model.  The  long-channel  strong-inversion 

transcapacitances  (see  Figure  4.26)  also  show  the  same  smooth  transition  across  the 
linear/saturation  boundary,  and  the  gate  and  bulk  transcapacitances  are  symmetric 
with  respect  to  the  source  and  drain  for  V^j  =  0  as  expected.  The  UFET  capacitance 

curves  show  good  qualitative  and  quantitative  agreement  with  the  corresponding 
results  from  MEDICI  presented  in  Figures  4.27  and  4.28. 

The  UFET  intrinsic  gate  capacitance  and  transcapacitances  for  the  short- 
channel  device  are  shown  in  Figures  4.29  and  4.30,  and  they  still  exhibit  smooth 
transitions  between  the  different  regions  of  operation.  The  model  correctly  predicts 
the  differences  in  the  capacitance  resulting  from  decrease  in  channel  length.  The 
simulation  results  of  Figure  4.29  show  the  effect  of  DIBL  in  weak  inversion,  and  they 
also  exhibit  the  expected  velocity  saturation  characteristics  in  strong  inversion  as 
opposed  to  the  pinch-off  behavior  of  a  long-channel  device.  As  in  the  long-channel 
case,  there  is  still  good  qualitative  agreement  between  the  UFET  capacitances  and 
the  MEDICI  results  shown  Figures  4.31  and  4.32.  The  quantitative  error  in  the  short- 
channel  devices  is  due  to  the  overlap  (parallel  plate)  and  fringing  (neglected) 
capacitance  modeling  in  UFET.  The  effect  of  DIBL  on  the  gate  capacitance  is  also 
over  predicted  which  is  due  to  the  assumption  of  the  parabolic  potential  along  the 
channel  to  define  the  two-dimensional  effects. 
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Figure  4.25  UFET-simulated  gate  capacitance  for  L^ff  =  9.95nm  with  Vps  =  OV. 
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Figure  4.26  UFET-simulated  transcapacitances  for  L^ff  =  9.95nm  with  Vgs  =  OV. 
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Figure  4.27  MEDICI-simulated  gate  capacitance  for  L^ff  =  9.95|im  with  Vgs  =  OV. 
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Figure  4.28  MEDICI-simulated  transcapacitances  for  L^ff  =  9.95|xm  with  Ygj  =  OV. 
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Figure  4.29  UFET-simulated  gate  capacitance  for  L^ff  =  0.3nm  with  Vgs  =  OV. 
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Figure  4.30  UFET-simulated  transcapacitances  for  L^ff  =  03\im  with  W^s  -  0^- 
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Figure  4.31  MEDICI-simulated  gate  capacitance  for  L^ff  =  0.3|im  with  Vgs  =  OV. 
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Figure  4.32  MEDICI-simulated  transcapacitances  for  L^ff  =  0.3nm  with  Vgs  =  OV. 


129 

Although  MEDICI  does  not  account  for  quantum  mechanical  effects  and 
polysilicon-gate  depletion  was  not  included  in  the  simulations,  the  impact  of  these 
two  effects  on  the  gate  capacitance  curves  can  be  shown  by  UFET.  In  Figure  4.33, 
the  three  curves  represent  the  gate  capacitance  without  polysilicon-gate  depletion  or 
quantum  mechanical  effects,  with  polysilicon-gate  depletion  and  with  polysilicon- 
gate  depletion  and  quantum  mechanical  effects.  Neither  of  these  effects  significantly 
influences  the  results  in  weak-inversion,  but  both  tend  to  increase  the  threshold 
voltage  and  reduce  the  capacitance  in  strong  inversion.  The  quantum  mechanical 
effects  in  the  accumulation  layer  also  lead  to  a  reduction  in  the  gate  capacitance 
below  flatband  similar  to  the  strong  inversion  results. 

4.7  Application  to  0.25|xm  Technology 

An  Semiconductor  Research  Corporation  (SRC)  member  company 
provided  measured  data  from  a  0.25|j.m  technology  (K.  Stiles,  private 
communication,  1997)  still  in  the  design  stage.  This  technology  has  a  5nm  gate  oxide 
with  an  LDD/LDS  structure  designed  for  a  2.5V  supply.  Since  the  process  was  still 
under  design  with  the  previous  generation  of  equipment,  the  process  control  of  AL  as 
a  function  of  channel  length  was  poor  (as  reported  by  the  company).  Some  liberty 
was  assumed  in  tuning  the  parameters  which  resulted  in  a  reasonable  fit  for  the  longer 
channel  devices  (as  shown  in  Figure  4.34),  but  the  physical  significance  of  the 
parameters  was  compromised  by  this  process.  The  significant  problems  seen  in  the 
subthreshold  plot  are  due  to  assuming  that  DL  was  a  constant,  and  the  preliminary 
accounting  for  RSCE/halo  effects  as  presented  in  Chapter  3.  Due  to  this  problem, 
effort  was  focused  solely  on  the  long-channel  device,  and  the  only  other  results 
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Figure  4.33  UFET-simulated  impact  of  polysilicon-gate  depletion  and  quantum 
mechanical  effects  on  gate  capacitance. 
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Figure  4.34  Subthreshold  plot  of  UFET-simulated  (lines)  and  measured  (points)  data 
in  0.25^m  teclinology. 
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presented  are  for  this  channel  length.  UFET  was  able  to  accurately  predict  the  both 
the  body  effect  (shown  in  Figure  4.35)  and  the  drain  bias  characteristics  off 
Figure  4.36.  As  in  the  simulated  0.35|im  technology,  the  source  of  the  error  in  the 
drain  bias  plot  is  due  to  the  constant  depletion  capacitance  assumption.  In  addition, 
the  company  also  supplied  data  taken  over  a  range  of  temperatures,  which  allowed 
the  evaluation  of  the  temperature  dependencies  in  the  model.  The  results  shown  in 
Figure  4.37  validate  the  temperature  modeling  in  UFET,  and  the  discrepancies  seen 
at  the  lower  temperatures  may  be  due  measurement  problems  (K.  Stiles,  private 
communication,  1997). 

4.8  Application  to  O.lS^m  Technology 
i 
4.8.1   0.15|im  Technology  Parameters 

Measured  data  was  supplied  by  an  SRC  member  company  (Y.  Taur, 

private  communication,  1997)  from  their  aggressive  0.15|im  retrograde  technology 

with  a  27A  gate  oxide,  LDD/LDS  structures  having  a  peak  doping  of  8x10     cm'^, 

a  peak  channel  doping  of  1x10  cm"^  and  halo  implants  processed  using  x-ray 
lithography.  Unfortunately,  the  measured  data  provided  shows  signs  of  pervasive 
contact  problems  and  device-to-device  variations.  Due  to  the  contact  resistance 
problem,  the  parameters  were  tuned  from  the  0.175|j.m  and  10|im  devices,  and  the 
0.15nm  device  was  modeled  by  varying  only  the  parasitic  source/drain  resistance. 
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Figure  4.35  Subthreshold  plot  of  UFET-simulated  (lines)  and  measured  (points)  data 
in  0.25|im  technology  showing  body  effect. 
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Figure  4.36  Linear-saturation  plot  of  UFET-simulated  (lines)  and  measured  (points) 
data  in  0.25|J.m  technology. 
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Figure  4.37  Thresholdplot  of  UFET-simulated  (lines)  and  measured  (points)  data  in 
0.25|im  technology  showing  temperature  effects. 
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4.8.2  Verification 

The  resulting  subthreshold  plots  are  shown  in  Figure  4.38  with  the 
subthreshold  slope  tuned  to  the  value  just  below  threshold.  Although  the  results  for 
the  long-channel  devices  are  reasonable,  the  short-channel  devices  show  significant 
discrepancies.  These  discrepancies  can  be  attributed  to  the  onset  of  punch-through 
current  and  to  the  simple  model  being  used  to  model  RSCE/halo  effects.  The 
threshold  plot  of  Figure  4.39  shows  that  UFET  is  able  to  model  the  characteristics  of 
the  measured  devices  in  a  deep-submicron  technology  for  the  long-channel  devices. 
However,  the  RSCE/halo  effects  cause  a  shift  in  the  threshold  voltage  which  is  not 
well  predicted  by  UFET,  indicating  that  the  RSCE/halo  model  in  UFET  does  not 
completely  account  for  this  effect.  (Note  that  the  errors  in  the  channel  current  for  the 
longer  channel  devices  are  within  the  10%  variation  seen  between  adjacent  devices 
with  the  same  W/L.)  The  drain-induced  effects  are  shown  in  Figure  4.40  show  an 
increased  error  in  subthreshold  slope  and  current  due  to  RSCE/halo  effects  and 
possible  punch-through  current  at  the  higher  drain  bias.  The  drain  bias  characteristics 
of  a  long-  and  short-channel  device  are  shown  in  Figures  4.41  and  4.42  respectively. 
Here,  the  long-channel  device  shows  significant  error  (-15%)  which  can  again  be 
attributed  to  the  model  assumption  of  a  constant  depletion  capacitance  along  the 
channel.  The  error  in  the  short-channel  device  can  be  attributed  to  velocity  overshoot 
(-10%  [Jan97])  which  is  not  modeled  and  to  possible  errors  in  tuning  the  value  of 
DL.  The  charge  modeling  has  also  been  further  verified  with  the  measured  gate 
capacitance  of  a  long-wide  device  as  shown  in  Figure  4.43.  Note  that  the  capacitance 
is  predicted  well  by  UFET  calibrated  only  via  DC  data. 
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Figure  4.38  Subthreshold  plot  of  UFET-simulated  (lines)  and  measured  (points)  data 
in  0. 1  S^im  technology. 
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Figure  4.39  Threshold  plot  of  UFET-simulated  (lines)  and  measured  (points)  data  in 

0.15^m  technology. 
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Figure  4.40  Subthreshold  plot  of  UFET-simulated  (lines)  and  measured  (points)  data 
in  O.lSlim  technology  showing  the  drain-induced  effects. 
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Figure  4.4 1   Linear-saturation  plotof  UFET-simulalcd  (lines)  and  measured  (points) 
data  in  0.15)xm  technology  for  L  =  lOjim. 
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Figure  4.42  Linear-saturation  plot  of  UFET-simulatcd  (lines)  and  measured  (points) 
data  in  0.15|lm  technology  for  L  =  0. 15|im. 
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Figure  4.43  Gate  Capacitance  plot  of  UFET-simulated  (line)  and  measured  (points) 
data  in  0.15|ini  technology  for  a  long-wide  device. 
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The  discrepancies  in  the  low-Vps  plots  in  this  0.15)lm  technology  are  due 
to  the  Reverse  Short-Channel  Effect  and  the  presence  of  a  halo.  The  preliminary 
accounting  in  UFET  for  these  effects,  which  are  important  in  scaled  technologies,  is 
being  done  using  a  model  first  developed  for  threshold  voltage.  This  accounting  for 
RSCE/halo  fails  to  model  this  dependence  properly,  and  this  is  an  area  for  further 
review  (see  Chapter  6).  The  discrepancies  in  the  drain  characteristics  seen  in  the 
long-channel  device  can  be  attributed  to  the  assumption  of  a  constant  depletion 
capacitance  along  the  channel.  Although  this  assumption  was  valid  for  the  two  other 
technologies  discussed,  it  breaks  down  for  this  0.15|im  technology.  This  indicates 
that  the  retrograde  doping  profile  may  not  be  as  steep  in  this  technology  which  allows 
the  depletion  depth  to  vary  more  along  the  channel  increasing  the  current  drive  in 
long-channel  devices  where  pinch-off  controls  the  characteristics.  The  need  for  an 
accounting  of  velocity  overshoot  also  needs  to  be  investigated  although  "clean"  data 
would  be  needed  for  this.  As  an  initial  assessment  of  this  effect,  hydrodynamic 
simulations  should  be  run  using  MEDICI  to  investigate  the  impact  of  velocity 
overshoot  in  extremely  scaled  devices  with  regard  to  the  UFET  formalism. 

4.9  Conclusion 

The  new  UFET  model  was  verified  against  MEDICI  simulations  and 
measured  data  for  the  DC  channel  current  across  all  biases  and  channel  lengths  with 
(overall)  good  results.  The  charge  modeling  was  also  verified  through  the  intrinsic 
capacitances  for  both  long-  and  short-channel  devices  across  all  biases  with 
quantitative  verification  of  the  gate  charge  model.  The  model  is  useful  not  only  for 
digital  design,  but  also  for  analog  design  where  the  derivatives  of  the  output 
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characteristics  are  important.  The  limited  number  of  physical,  process-based  model 
parameters  with  inherent  correlation  between  the  parameter  values  and  device 
characteristics  also  make  the  model  ideally  suited  for  process  design/integration  and 
statistical  simulations. 


CHAPTER  5 
MEDICI  BENCHMARK  AND  A  SENSITIVITY  STUDY 


5.1   Introduction 

Two  important  issues  for  any  device  model  are  parameter  extraction/ 
evaluation  and  correlation  to  process  parameters  as  needed  for  process  design/ 
integration,  statistical  simulations  and  sensitivity  analyses.  In  models  such  as 
BSIM3v3  with  a  large  number  of  parameters,  a  copious  amount  of  data  across  several 
channel  lengths  and  all  terminal  biases  is  needed  in  order  to  extract  the  parameters 
controlling  the  length  and  bias  dependencies.  As  a  result  of  the  excessive 
optimization  needed  for  such  models,  the  correlation  between  model  parameters  and 
process  variations  may  not  be  clear,  may  be  spread  across  several  parameters  or  may 
not  exist  at  all.  With  its  small  number  of  physical  parameters,  the  process-based 
UFET  model  makes  parameter  evaluation  relatively  easy  as  demonstrated  in 
Chapter  4  while  maintaining  a  strong  correlation  to  process  variations. 

In  this  chapter,  UFET  is  first  applied  to  a  simulated  0.1  |xm  technology  to 
demonstrate  its  applicability  to  extremely  scaled  technologies.  Based  on  only  DC 
calibration  involving  minimal  data,  transient  and  AC  figures  of  merit  are  simulated, 
and  their  sensitivities  are  examined.  A  circuit  application  is  then  presented  using  the 
0.15|J.m  technology  in  Chapter  4  to  demonstrate  the  UFET's  utility  in  circuit 
sensitivity  analyses. 


145 


146 

5.2  MEDICI  Benchmark 

In  order  to  further  demonstrate  the  utility  of  this  new  process-based 
model,  a  O.lum  technology  was  simulated  using  MEDICI,  and  the  UFET  parameters 
were  tuned  using  a  limited  set  of  DC  data  from  two,  short-  and  long-channel  (L  = 
0.  l^m  and  lOjlm),  devices  over  a  limited  range  of  terminal  biases: 

•  Ich-Vgs  for  Vds  =  0. 1 V  and  Vdd 

•  IcH-^DS  for  Vgs  =  0-5V,  1  .OV  and  1 .5V 

•  IcH-Vosfor  Vi3s  =  0.1Vand  Vbs  =  OV  and -1.5V 

Trying  to  use  a  model  such  as  BSIM3v3  for  this  exercise  is  problematic,  if  not 
impossible,  due  to  the  copious  amounts  of  data  required  for  the  intensive 
optimization  required  for  parameter-extraction  (K.  R.  Green,  private  communication, 
1997).  After  tuning  the  model  parameters,  UFET-  and  MEDICI-simulated  results 
were  compared  for  two  intermediate  channel  lengths  (L  =  0.3|im  and  l|im)  over  the 
same  range  of  terminal  biases.  In  addition,  transient  and  AC  simulations  were  done 
with  MEDICI  and  UFET  to  find  the  intrinsic  propagation  delay  (tp^)  and  cutoff 
frequency  (f-j-)  of  the  four  devices.  In  order  to  show  the  physical  nature  of  the  model, 
sensitivity  simulations  were  performed  with  assumed  variations  in  t^^  and  AL  from 
the  nominal  values  to  assess  the  impact  on  device  performance. 

5.2.1  O.lum  Technology 

The  technology  was  designed  to  achieve  a  0.5V  threshold  voltage  with 
minimal  threshold  voltage  roll-off  due  to  charge-sharing  and  drain  induced  effects. 
Based  on  published  0.1|j,m  technologies  [Iwa93,  Mii94],  the  gate  oxide  was  chosen 
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17  -3 

to  be  5nm  with  a  substrate  doping  of  7x10  cm  ,  a  junction  depth  of  0.08|Xm 
corresponding  to  the  extrinsic  source/drain  junction  depth  and  a  gate  overlap  of 
0.01 5|xm  due  to  the  lateral  diffusion  of  the  source  and  drain  junctions.  The  peak 

18  -3 

doping  was  set  at  2x10  cm  with  a  Gaussian  profile  near  the  surface  determined 
by  the  threshold  voltage  criterion.  In  order  to  reduce  simulation  time  several 
simplifications/assumptions  were  made  about  the  device  structure: 

•  The  contact  resistance  was  assumed  to  be  zero. 

•  The  polysilicon-gate  doping  was  assumed  to  be  large, 
i.e.,  no  polysilicon-gate  depletion  effect. 

•  The  lateral  doping  profile  was  assumed  to  be  uniform, 
i.e.,  no  RSCE. 

•  The  simulation  was  done  using  only  drift-diffusion  variables, 
i.e.,  no  velocity  overshoot. 

Since  the  device  simulations  were  done  using  an  older  version  of  MEDICI,  the 

quantum  mechanical  effects  were  not  included  in  this  benchmark.  A  summary  of  the 

technology  and  simulation  parameters  is  shown  in  Table  5.1. 

5.2.2  DC  Results 

After  the  DC  MEDICI  simulations  were  completed,  the  model  parameters 
were  tuned  with  the  0.1  flm  and  10|j.m  devices  using  the  procedure  outlined  in 
Chapter  4,  and  this  set  of  model  parameters  was  then  used  to  predict  the  0.3|im  and 
Ijim  device  characteristics.  Note  that  although  the  model  parameters  were  evaluated 
using  only  the  0.1  |im  and  10)xm  devices,  the  results  are  presented  using  all  four 
channel  lengths  in  some  of  the  figures  to  follow. 
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Table  5.1   Nominal  0.1  |im  Technology 


Parameter 

Value 

tox 

5nm 

Xj 

0.08nm 

AL 

0.03|lm 

Substrate  Doping 

7xl0"   cm"' 

Peak  Doping 

2xl0'*  cm"' 

Vdd 

1.5V 

V,(Vds  =  0.1V) 

0.5V 

Reverse  Short-Channel  Effect 

No 

Polysilicon-gate  Depletion 

No 

LDD 

No 

Contact  Resistance 

No 

Hydrodynamic  Simulation 

No 

Quantum  Mechanical  Effect 

No 
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In  the  low-Vps  subthreshold  plot  of  Figure  5.1,  UFET  accurately  models 
the  subthreshold  slope  and  current  of  all  four  channel  lengths  with  no  substrate  bias. 
UFET  also  models  the  behavior  above  threshold  as  shown  in  Figure  5.2  again  with 


channel  lengths  are  presented  in  Figure  5.3.  The  long-channel  device  shows  no  drain- 
induced  effects  as  expected,  and  the  DIBL  effects  in  the  subthreshold  region  of 
operation  increase  as  the  channel  length  becomes  smaller.  There  is  good  agreement 
between  the  UFET-  and  MEDlCI-predicted  results  with  no  tuning  of  the  DIBL  effect 
for  the  intermediate  channel  lengths.  The  DIBL  effect  in  the  0.3fxm  device  are  over- 
predicted  which  may  mean  that  the  drain-induced  effects  tend  to  diminish  faster  than 

the  —  dependence  in  the  model  derived  by  assuming  a  parabolic  potential  for  the 

analysis  of  the  two-dimensional  effects. 

The  body  effect  is  shown  in  Figures  5.4-5.7,  and  there  is  good  agreement 
for  the  three  longest  channel  lengths  between  MEDICI  and  UFET.  However,  the 
MEDICI  results  for  the  0.1  (tm  device  show  an  error  due  to  the  severe  two- 
dimensional  effects  in  this  device  due  to  the  overlapping  source  and  drain  depletion 
regions  (which  was  later  verified  using  MEDICI).  This  implies  that  the  technology 
is  not  optimally  designed  since  punch-through  would  occur  at  higher  drain  biases, 
but  it  is  important  to  note  that  UFET  was  still  able  to  predict  the  short-channel  effects 
until  the  overlapping  depletion  regions  completely  dominated  the  characteristics. 

The  drain  characteristics  are  shown  in  Figures  5.8-5.1 1  with  good 
agreement  (less  than  10%  error  for  the  intermediate  channel  lengths)  between  UFET 
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Figure  5. 1   Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
in  nominal  0.1  |im  technology. 
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Figure  5.2  Threshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data  in 
nominal  0.1|a.m  technology. 
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Figure  5.3  Subthreshold  plots  of  UFET-(lines)  and  MEDICl-(points)simulated  data 
in  nominal  0.1)xm  technology  showing  DIBL. 
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Figure  5.4  Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
in  nominal  0.1|im  technology  for  L  =  0.1|im  showing  V3S  dependence. 
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Figure  5.5  Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
in  nominal  0.1|xm  technology  for  L  =  0.3nm  showing  Ygs  dependence. 
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Figure  5.6  Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
in  nominal  0.1)i,m  technology  for  L  =  l|lm  showing  Ygj  dependence. 
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Figure  5.7  Subthreshold  plots  of  UFET-(lines)  and  MEDICI-(points)simulated  data 
in  nominal  0.1  |im  technology  for  L  =  10|lm  showing  Vgj  dependence. 


157 


'-aa.e  =  0-''^im     \N  =  ^O^lm 


-gate 


•  v„ 

=  1.5V 

■  v,. 

=  1.25V 

♦Vos 

=  t.OV 

*Vos 

=  0,75V 

TV,, 

=  0.5V 

▼  V„ 

=  0.25V 

0     »*♦?    TTTTTTT^ 

0.0  0.5 


1.0 


T    T    T    T 

1.5 


Vos      (V) 


Figure  5.8  Triode-Saturation  plots  of  UFET-(lines)  and  MEDICI-(points)simulated 
data  in  nominal  0.1  urn  technology  for  L  =  O.I|lm. 
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Figure  5.9  Triode-Saturation  plots  of  UFET-(lines)  and  MEDICI-{points)simulated 
data  in  nominal  0.1  nm  technology  for  L  =  0.3|im. 
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Figure  5.10  Triode-Saturation  plots  of  UFET-(lines)  and  MEDICI-(points)simulated 
data  in  nominal  0.1)ini  technology  for  L  =  l|im. 
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Figure  5.11  Triode-Saturation  plots  of  UFET-(lines)  and  MEDICI-(points)simulated 
data  in  nominal  0.1  urn  technology  for  L  =  lOnm. 
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and  MEDICI.  There  is  a  good  match  for  the  channel  current  with  only  a  small  error 
in  the  drain  conductance  for  the  0.1  |im  device.  As  was  the  case  for  the  0.35nm 
technology,  the  likely  cause  for  the  error  in  the  drain  conductance  is  the  DL  assumed. 
Since  the  technology  was  simulated  with  MEDICI,  the  DL  was  assumed  to  represent 
the  difference  between  the  physical  gate  length  and  the  metallurgical  channel  length, 
but  a  recent  paper  [Niu97]  indicates  that  the  effective  (i.e.,  gate-controlled)  channel 
length  is  larger  than  the  metallurgical  channel  length,  suggesting  that  DL  should  be 
smaller  than  the  value  used.  This  will  decrease  the  drain  conductance  while 
necessitating  an  increase  in  VSAT  which  tends  to  be  different  from  the  MEDICI 
value  due  to  the  differing  mobility  models  used.  The  increase  in  VSAT  will  tend  to 
decrease  the  error  seen  in  the  0.3|lm  device.  For  the  l[im  and  10p,m  devices,  the 
dominant  source  of  error  is  the  constant  depletion  capacitance  approximation  used. 


5.2.3  Transient  and  AC  Results 

In  order  to  demonstrate  the  predictive  nature  of  UFET  in  a  digital 
simulation,  the  gate  delay  of  the  devices  was  investigated.  Due  to  limitations  in  the 
simulation  capabilities  of  generic  MEDICI  (i.e.,  without  the  Circuit  Analysis 
Advanced  Application  Module),  only  the  intrinsic  pull-down  delay  of  the  devices 
was  simulated  as  shown  in  Figure  5.12.  Also  due  to  the  poor  time-step  control  which 
required  the  use  of  a  small  constant  time-step  in  these  simulations,  the  intrinsic  delay 
of  the  10|im  devices  was  not  simulated.  As  shown  in  Figure  5.13,  UFET  accurately 
predicts  the  channel  length  and  supply  voltage  dependencies  of  the  intrinsic  delay  for 
the  three  different  channel  lengths.  The  slight  differences  between  the  MEDICI  and 
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Figure  5.12  Intrinsic  delay  circuit. 
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Figure  5.13  UFET-(lines)  and  MEDICI-(points)simulated  intrinsic  pull-down  delay 
in  nominal  0.1  |im  technology. 
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UFET  results  can  be  explained  by  the  differences  in  saturated  channel  currents  seen 
in  Figures  5.8-5.1 1. 

Application  of  UFET  for  analog  design  requires  verification  of  the 
(derivatives  of  the)  channel  current  and  terminal  charges.  This  has  been  done  by 
simulating  the  cutoff  frequency  [Gra93], 


with  both  MEDICI  and  UFET.  In  Figure  5.14,  the  cutoff  frequency  is  shown  as  a 

function  of  both  channel  length  and  gate  bias,  and  UFET  accurately  predicts  both 

trends.  The  maximum  cutoff  frequency,  which  occurs  at  approximately  gni(max)  ^^ 

can  be  seen  in  (5.1),  can  be  extracted  from  the  data  in  Figure  5.14  and  plotted  as  a 

function  of  channel  length  as  shown  in  Figure  5.15.  UFET  is  able  to  accurately  model 

the  cutoff  frequency  over  two  orders  of  magnitude  in  channel  length. 

The  results  for  the  intrinsic  delay  and  cutoff  frequency  were  obtained 

using    the    parameter    set    calibrated    from    limited    DC    data.    This    shows    the 

predictability  of  UFET,  which  has  practical  applications  in  assessing  preliminary 

data  from  new  technologies.  The  physical  nature  of  the  model  parameters  also  allows 

UFET  to  be  used  in  investigations  of  device  sensitivity  to  process  variations  as  will 

now  be  demonstrated. 

5  '•■' 
5.2.4  Device  Sensitivity 

As  both  vertical  and  lateral  device  dimensions  continue  to  shrink,  it  is 

becoming  even  more  important  to  tightly  control  the  fabrication  process  in  order  to 
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Figure  5.14  UFET-(lines)  and  MEDICI-(points)simulated  cutoff  frequency  in 
nominal  0.1  )xm  technology. 
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Figure  5.15  UFET-(line)  and  MEDICI-(points)simulated  cutoff  frequency  in 
nominal  0.1  [xm  technology. 


167 

minimize  device  variations  across  lots  (to  improve  yield)  and  also  wafers  (in  order 
to  ensure  the  device  matching  needed  in  some  circuits).  At  the  device  level, 
sensitivity  to  variations  in  oxide  thickness  and  lateral  diffusion  (AL)  can  be  easily 
shown  using  the  process-based  model  parameters  in  UFET. 

The  oxide  thickness  was  assumed  to  have  a  maximum  10%  variation 
which  is  consistent  with  published  variations  [Aro95].  For  the  case  of  pull-down 
delay,  the  variation  in  oxide  thickness  is  a  first-order  effect  due  to  the  linear  variation 
in  drive  current  with  to^.  As  shown  in  Figures  5.16  and  5.17,  UFET  models  the 
increase  and  decrease  in  the  pull-down  delay  associated  with  the  change  in  tg^.  The 
cutoff  frequency  as  given  by  (5.1)  is  independent  of  t^^  to  first-order,  and  the 
resulting  plots  of  fj  in  Figure  5.18  show  the  device  has  little  sensitivity  to  variations 
in  the  oxide  thickness. 

Assuming  a  typical  variation  in  AL  of  50%  [Aro95],  sensitivity 
simulations  of  l„^  and  fj  were  performed.  As  expected  from  a  simple  drive  current 
analysis,  the  variation  of  the  pull-down  delay  is  greatest  for  short-channel  devices 
with  almost  no  sensitivity  in  the  long-channel  case  as  shown  in  Figures  5.19  and 
5.20.  A  similar  result  for  the  cutoff  frequency  is  shown  in  Figure  5.21  with  the  UFET 
model  predicting  the  trends  well. 


5.3  Operational  Amplifier  Sensitivity  Study 

In  order  to  demonstrate  the  utility  of  UFET  in  a  circuit  application,  an 
operational  amplifier  [A1187]  is  examined  to  determine  the  variation  in  circuit 
performance  in  relation  to  process  variations.  The  CMOS  operational  amplifier 
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Figure  5.16  UFET-(lines)  and  MEDICI-(points)simulated  intrinsic  pull-down  delay 
in  O.lum  technology  with  tg^  =  5.5nm. 
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Figure  5.17  UFET-(lines)  and  MEDICI-(points)simulated  intrinsic  pull-down  delay 
in  O.lum  technology  with  t^^  =  4.5nm. 
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Figure  5.18  UFET-(lines)  and  MEDICI-(points)simulated  cutoff  frequency  in  0. 1  |im 
technology  with  10%  oxide  thickness  variations. 
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Figure  5.19  UFET-(lines)  and  MEDICI-(points)simulated  intrinsic  pull-down  delay 
in  0.1  nm  technology  with  AL  =  0.045|lm. 
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Figure  5.20  UFET-(lines)  and  MEDICI-(points)simulated  intrinsic  pull-down  delay 
in  0.1|Xm  technology  with  AL  =  0.015|im. 
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Figure  5.21  UFET-(lines)  and  MEDICI-{points)simulated  cutoff  frequency  in  0. 1  ^im 
technology  with  AL  50%  variations. 
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designed  using  normalized  channel  lengths  and  presented  at  a  SEMATECH 
workshop  (B.  M.  Tenbroek,  private  communication,  1997)  is  chosen;  the  schematic 
is  shown  in  Figure  5.22.  Since  UFET  is  process-based,  the  sensitivity  study  is  a 
relatively  straightforward  process  since  the  model  parameters  are  strongly  correlated 
with  the  process  variables. 

The  operational  amplifier  is  a  two-stage  folded  cascode  design.  The  first 
(input)  stage  consists  of  a  differential  PMOS  pair  with  an  active  load.  The  cascode 
circuit  increases  the  gain  of  the  first  stage  by  increasing  the  output  resistance  seen  by 
the  first  stage.  The  second  stage  (output)  stage  is  a  simple  class-A  amplifier.  The 
compensation  is  provided  by  the  capacitor  in  series  with  a  PMOS  transistor  which 
cancels  the  non-dominant  pole  to  improve  the  phase  margin  of  the  circuit  [A1187].  In 
addition,  there  are  several  NMOS  transistors  which  function  as  the  bias  circuit  for 
the  operational  amplifier.  I 


5.3.1   Process  Variations 

Using  the  0.15|im  technology  presented  in  Chapter  4,  a  subset  of  the 
process  variables  (to^,  ALnmos-  ^LpMos-  Npqly.nmos  and  NpoLY.PMOs)  is  selected. 
Of  these,  the  largest  uncontrolled  process  variation  is  the  lateral  diffusion  which 
affects  the  threshold  voltage  in  extremely  scaled  devices  and  the  overlap 
capacitances  which  become  more  important  in  extremely  scaled  devices.  Variations 
in  oxide  thickness  and  polysilicon-gate  doping  concentration  also  lead  to  variations 
in  threshold  voltage.  I  : 

Although  only  NMOS  devices  were  calibrated,  the  PMOS  devices  are 
assumed  to  be  "symmetric"  with  a  60%  reduction  in  the  mobility,  i.e.,  all  model 
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parameters  are  the  same  except  for  a  mobility  reduction.  In  addition,  the  contact 
resistance  problem  seen  in  the  data  is  ignored,  and  the  same  model  parameters  are 
used  for  all  channel  lengths.  Based  on  published  process  variations  [Rio95,  Aro95], 
variations  in  the  process  variables  of  interest  are  set.  The  oxide  thickness  is  presumed 
to  vary  by  10%  from  the  nominal  process  value  of  2.7nm,  and  the  lateral  diffusion 
lengths  for  both  devices  are  assumed  to  vary  by  50%  from  the  nominal  value  of 
0.03|J.m.  The  influence  of  the  polysilicon-gate  depletion  effect  on  the  circuit 
performance  is  examined  by  varying  the  polysilicon-gate  doping  concentrations  from 

1x10  cm"  (severe  gate  depletion  effects)  to  1x10  cm  (negligible  gate  depletion 
effects). 

5.3.2  Sensitivity  Results 

Using  the  presumed  process  variations,  three  figures  of  merit  important  to 
analog  circuit  designers  (DC  gain,  unity  gain  frequency  and  phase  margin)  are 
examined  to  determine  their  dependencies  on  process  variations.  The  results  could 
then  be  used  by  a  circuit  designer  to  determine  if  the  operational  amplifier  needs  to 
be  redesigned  in  order  to  better  tolerate  the  process  variations. 

Using  the  nominal  technology  channel  (L  =  0.15|xm)  with  the  nominal 
process  values  results  in  an  operational  amplifier  with  a  good  DC  gain  (75dB),  unity 
gain  frequency  (42MHz)  and  phase  margin  (90°)  as  shown  by  the  UFET  simulation 
results  in  Figures  5.23  and  5.24.  Although  these  results  indicate  that  the  operational 
amplifier  may  be  over-compensated  as  implied  by  the  large  phase  margin,  the  pole- 
zero  cancellation  circuit  is  not  optimally  designed  since  the  channel  resistance  is 
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Figure  5.23  UFET-simulated  open  loop  gain  for  nominal  0.15|lm  technology. 
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Figure  5.24  UFET-simulated  open  loop  phase  for  nominal  0.15jlm  technology. 
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smaller  than  anticipated  due  to  the  effects  of  AL.  Examination  of  the  correlation 
between  phase  margin  and  the  process  variations  may  indicate  a  reduction  in  the 
value  of  the  compensation  capacitor  is  possible  which  will  reduce  the  circuit  area  and 
increase  the  slew  rate  [A1187]  after  the  circuit  has  been  redesigned. 

Over  the  process  variations  in  the  gate  oxide  thickness,  the  DC  gain  shows 
a  small  sensitivity  as  illustrated  in  Figure  5.25  with  a  larger  sensitivity  for  the  unity 
gain  frequency.  Since  the  main  purpose  of  the  NMOS  transistors  is  to  control  the  bias 
conditions  (as  forced  by  Ibj^s  in  Figure  5.22),  there  is  little  sensitivity  to  the  process 
variations  associated  with  the  NMOS  devices.  This  is  shown  for  variations  of  NMOS 
lateral  diffusion  in  Figure  5.26  and  for  NpoLv  '"  Figure  5.27.  The  PMOS  transistors 
control  the  dominant  portions  of  the  operational  amplifier  circuit,  and  it  is  expected 
that  there  will  be  a  larger  sensitivity  for  these  devices  than  for  the  NMOS.  This  is 
true  for  both  the  PMOS  lateral  diffusion  and  PMOS  polysilicon-gate  doping  as 
shown  in  Figures  5.28  and  5.29. 

Using  the  results  of  this  sensitivity  study  allows  a  series  of 
recommendations  to  be  made  regarding  design  improvements  for  this  circuit.  The 
circuit  is  relatively  insensitive  to  t^^  and  AL^mos  variations  in  the  process.  It  is  also 
relatively  insensitive  to  NpoLY.NMOS  ^"'^  '^POLY.PMOS  ^°^  doping  concentrations 
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larger  than  3x10  cm  ,  which  are  easily  achieved  [Aro95].  However,  the  design  is 
sensitive  to  variations  in  the  PMOS  lateral  diffusion,  and  the  circuit  should  be 
redesigned  to  reduce  this  sensitivity.  This  is  accomplished  by  increasing  the  channel 
lengths  and  widths  of  the  PMOS  devices  controlling  the  circuit  performance.  This 
will  decrease  the  sensitivity  of  the  circuit  at  the  penalty  of  increased  area.  However, 


180 


a- 


-5 


-5  0 

t^^  Variation  (%) 


10 


Figure  5.25  UFET-simulated  t^^  sensitivity. 
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Figure  5.26  UFET-simulated  ALn^jos  sensitivity. 
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Figure  5.27  UFET-simulated  NpQLy.NMOS  sensitivity. 
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Figure  5.28  UFET-simulated  ALp^jos  sensitivity. 
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the  results  also  show  that  the  circuit  is  over-compensated,  and  the  increased  area  can 
be  partially  offset  by  reducing  the  size  of  the  compensation  capacitor.  The  resulting 
operational  amplifier  using  0.45nm  for  the  normalized  channel  lengths  in  the  PMOS 
devices  has  improved  DC  gain  (86dB),  a  smaller  unity  gain  frequency  (37MHz)  and 
a  similar  phase  margin  (87°),  but  more  importantly  the  circuit  sensitivity  to  ALp^os 
has  been  greatly  reduced  as  shown  in  Figure  5.30.  The  problem  with  the  pole-zero 
cancellation  has  also  been  eliminated  as  shown  in  the  gain  and  phase  plots 
Figures  5.31  and  5.32. 


5.4  Conclusion 

UFET's  small  number  of  process-based  parameters  simplify  the  parameter 
extraction/evaluation  procedure,  and  the  correlation  to  process  parameters  enable 
UFET  to  be  easily  used  in  modeling  device  sensitivity.  This  was  demonstrated  with 
a  MEDICI  exercise  where  UFET  was  used  to  simulate  an  extremely  scaled 
technology,  and  to  also  examine  device  sensitivity  to  variations  in  the  oxide 
thickness  and  lateral  diffusion  in  this  technology.  This  demonstration  was  then 
extended  by  using  UFET  to  examine  the  sensitivity  of  a  CMOS  operational  amplifier 
to  normal  process  variations.  The  results  show  that  UFET  can  be  a  viable  alternative 
to  the  continued  use  of  empirical  compact  models  for  both  circuit  simulation  and 
device  sensitivity  studies. 


186 


-50    -40    -30    -20-10      0       10     20     30     40     50 

^LpMos  Variation  (%) 


Figure  5.30  UFET-simulated  ALp^jos  sensitivity  in  the  redesigned  circuit. 
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Figure  5.31   UFET-simulated  open  loop  gain  for  redesigned  0.1 5nm  technology. 
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Figure  5.32  UFET-simulated  open  loop  phase  for  redesigned  0.15|xm  technology. 


CHAPTER  6 
CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 


6.1  Conclusions 

The  many  deficiencies  of  existing  compact  short-channel  MOSFET 
models  imply  the  need  for  a  new  physical,  predictive  model.  This  dissertation 
presented  such  a  model  for  deep-submicron  devices  (UFET)  with  C„  continuity  that 
relies  on  approximate  two-dimensional  analyses,  which  are  computationally 
efficient,  to  account  for  the  predominant  short-channel  effects.  These  solutions  of 
two-dimensional  analyses  are  physically  linked  and  are  expressed  in  simple  terms 
that  yield  physical  insight  into  device  behavior  as  well  as  a  predictive  circuit  model 
with  utility  for  TCAD,  e.g.,  for  statistical  simulations  and  circuit  sensitivity  analyses, 
as  well  as  design.  The  parameter-evaluation  procedure  was  outlined,  and  the 
advantages  of  a  physical  model  where  the  process-based  model  parameters  are 
closely  tied  to  the  device  physics  and  structure  were  demonstrated.  This  is  in  contrast 
to  models  such  as  BSIM3  [Che96]  which  are  dependent  on  extensive  parameter 
optimization  that  can  lead  to  non-physical  parameter  sets.  Validation  of  UFET  was 
then  presented  using  device  simulations  and  measured  data,  and  the  utility  of  the 
model  was  shown  through  a  MEDICI  benchmark  and  a  sensitivity  analysis  of  a 
CMOS  operational  amplifier.  The  UFET  model  implemented  in  the  Type-I  interface 
is  available  from  the  UF  SOI  Group  at  the  University  of  Florida. 
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6.2  Recommendations  for  Future  Work 

Although  much  progress  was  made  in  the  development  of  the  UFET 
MOSFET  model,  there  is  still  much  more  work  to  be  done  in  order  to  make  it  the 
model  of  choice  for  modeling  deep-submicron  MOSFETs.  The  focus  of  this  work 
was  the  development  of  a  process-based  MOSFET  model  with  C^  continuity.  The 

next  phase  of  this  research  should  focus  on  the  effects  which  were  not  included  or 
not  well  modeled  in  UFET.  Along  with  upgrading  the  extrinsic  portions  of  the  model, 
this  will  make  UFET  a  complete  MOSFET  model  for  deep-submicron  technologies 
which  can  be  accepted  as  the  industry  standard. 

The  areas  for  further  UFET  investigation  are  listed  below  in  order  of 
importance 


RSCE/halo  effects  1 

Overlap  and  fringing  capacitances 

Narrow  width  effects 

Non-local  effects:  velocity  overshoot  and  impact  ionization 

Explicit  accounting  for  depletion  charge  variation  along  channel 

Leakage  current  and  junction  perimeters 

Tunneling  currents 

Noise  modeling  i 


The  big  hurdle  for  physical  MOSFET  modeling  is  the  Reverse  Short-Channel  Effect 
due  to  Transient  Enhanced  Diffusion  and  the  introduction  of  nonuniform  lateral 
doping  profiles,  i.e.,  halos,  which  result  in  a  changing  device  structure  as  a  function 
of  channel  length.  Although  a  preliminary  model  for  this  effect  was  implemented  in 
the  UFET  formalism,  the  comparisons  with  measured  data  suggest  that  the  RSCE 
model  should  be  improved. 

At  this  time,  all  of  the  effects  in  UFET  are  based  on  local-field  models. 
The  non-local  impact-ionization  model  in  UFSOI  [Fos97]  should  be  made  consistent 
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with  the  UFET  modeling  approach  in  order  to  account  for  the  non-local  effects  in 
extremely  scaled  devices. 

It  will  also  be  necessary  to  improve  the  modeling  of  the  extrinsic  portions 
of  the  device.  The  overlap  capacitances  are  modeled  as  simple  parallel-plate 
capacitors,  and  the  fringing  capacitances  are  ignored.  The  overlap  capacitances  are 
bias-dependent  and  should  be  modeled  as  such.  Due  to  significant  errors  in  the 
reverse-bias  junction  leakage  current,  an  approach  similar  to  that  in  the  UFSOI 
models  should  be  implemented  where  physical  models  for  generation  and 
recombination  are  used.  This  upgrade  should  also  account  for  the  perimeter 
components  in  the  leakage  current  and  the  junction  depletion  charges.  It  is  also 
necessary  to  account  for  the  narrow-width  effects  that  are  seen  in  measured  devices 
due  to  the  three-dimensional  effects  associated  with  the  isolation  structures.  The 
tunneling  currents  associated  with  the  scaled  gate  oxides  should  be  accounted  for, 
and  the  GIDL  model  currently  being  implemented  in  UFSOI  should  be  added  to 
UFET  when  the  development  is  complete. 

The  channel-current  analysis  can  also  be  upgraded  through  the  explicit 
integration  of  the  depletion  charge  along  the  channel  rather  than  assuming  a  constant 
depletion  capacitance.  Although  this  will  (slightly)  complicate  the  analysis  for 
Av|/|^',  the  trade-off  for  increased  accuracy  seems  reasonable. 

At  present,  UFET  lacks  a  noise  model.  The  thermal  noise  should  be 
modeled  as  proposed  by  Tsividis  [Tsi87]; 


Ai2 
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It  is  also  recommended  that  the  implementation  of  non-quasi-static  effects  be 
examined  in  order  to  confidently  use  the  models  for  high-frequency  analog  design. 
Although  the  UFET  model  has  been  moved  from  MMSPICE  into  the  API, 
where  the  model  is  not  completely  isolated  from  the  simulator,  this  process  needs  to 
be  continued  with  UFET  moving  along  to  the  next  interface  standard  (either  Type-II, 
Verilog,  MAST  or  some  similar  approach),  which  will  allow  the  model  to  be 
designed  without  regard  to  the  simulation  platform  used.  It  is  also  a  necessity  that  the 
parameter-evaluation  procedure  outlined  in  Chapter  4  be  implemented  in  software 
(e.g.,  UTMOST,  or  a  similar  package).  This  will  enable  UFET  to  be  more  easily 
distributed  and  used  effectively  and  pervasively  in  industry. 
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